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Abstract. In this paper wc propose a process of lagrangian reduction and 
reconstruction for nonholonomic discrete mechanical systems where the action 
of a continuous symmetry group makes the configuration space a principal 
bundle. The result of the reduction process is a discrete dynamical system 
that we call the discrete reduced system. We illustrate the techniques by 
analyzing two types of discrete symmetric systems where it is possible to go 
further and obtain (forced) discrete mechanical systems that determine the 
dynamics of the discrete reduced system. 

1. Introduction 

The elimination of degrees of freedom of a symmetric mechanical system, the ba- 
sic goal of reduction theory, is an old subject that dates back to the mid-nineteen 
century. The work of Routh in the context of abelian symmetries of classical me- 
chanical systems was extended by many others in an effort to explore different 
aspects of the reduction process. Currently, there arc well developed theories of 
reduction in the Hamiltonian setting, where the emphasis is in the reduction of the 
symplectic and Poisson structures, and the lagrangian setting, mostly focused on 
the reduction of variational principles. The literature in this area is vast; modern 
references are, for instance, [18, 2, 22, 8]. 

Discrete time mechanical systems have been considered in the literature since 
the 1960s, mostly as a way to approximate and model the behavior of (continuous) 
mechanical systems (see [21] and the references therein). The dynamics of nonholo- 
nomic discrete mechanical systems has been introduced by J. Cortes and S. Mar- 
tinez in [10]. Many characteristics of mechanical systems have a discrete analogue. 
One important feature whose discrete analogue has been exposed only partially in 
the literature is the reduction of symmetries. The purpose of the present work 
is to describe a reduction and reconstruction process for discrete time mechanical 
systems with nonholonomic constraints, in the lagrangian setting. 

There are many reasons for being interested in the reduction of a symmetric 
system, be it continuous or discrete. On the one hand, the reduced system has less 
degrees of freedom, which may be an advantage when trying to solve the equations 
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of motion. On the other, in many cases, the reduced system encodes the core 
dynamics of the system of interest. A good example of the advantage of working 
on the reduced system in the discrete case is presented by S. Jalnapurkar et aL 
in [IG] where they show how the reduced system is free from geometric phases' 
effects allowing them to observe interesting dynamical structures, which were hard 
to separate in the unreduced system. 

There are already several results on the reduction and reconstruction of sym- 
metric nonholonomic discrete mechanical systems available in the literature. The 
case where the configuration space of the system Q is a Lie group has been studied 
by Y. Fedorov and D. Zenkov in [13, 12] and by R. McLachlan and M. Perlmut- 
ter in [2.3]. They obtain reduced equations of motion on the Lie algebra of the 
Lie group of symmetries which coincide, in the unconstrained case, with those ob- 
tained by A. Bobcnko and Y. Suris in [4]. The case of Chaplygin systems, that 
is, when the symmetry group acts in such a way that no symmetry direction is 
compatible with the given constraints is analyzed by J. Cortes in [ 1 1] and, in con- 
text of groupoids, by D. Iglesias ct al. in [15]. Also, the unconstrained case has 
been treated by S. Jalnapurkar et al. in [ ] and, using discrete connections, by M. 
Leok et al. in [19]. A different approach using groupoids is given in [15]; in this 
case, the abstract reduction theory of groupoids produces a general reduced system 
on a groupoid; our approach differs from theirs, ours being more elementary — we 
stay in the "pair groupoid" case — and, at the same time, more explicit because we 
incorporate additional information (connections) to our models. 

The reduction results that we present here are modeled on the results of H. Ccn- 
dra et al. in [S] for continuous systems. They consider the action of a Lie group G 
on a configuration manifold Q so that the quotient map n : Q ^ Q/G is a principal 
bundle and the other data is G-invariant. They relate the variational principle that 
determines the dynamics of the original system with a reduced variational principle, 
that determines the dynamics of the reduced system. Eventually, equations of mo- 
tion are derived from both variational principles. The introduction of a connection 
on the principal bundle tt serves them to, first, construct an isomorphic model for 
the natural reduced space TQ/G and, also, to split the reduced variational principle 
and equations of motion into horizontal and vertical parts. Appropriate choice of 
a connection can lead to a simplified analysis of a specific mechanical system. 

At a philosophical level, a common approach to discrete mechanics consists of 
replacing the tangent bundle TQ hy Q x Q, with the idea that infinitesimal dis- 
placement — velocities — are replaced by finite displacement — pairs of points. Even 
though this is a powerful idea, there are some very important differences between 
TQ and QxQ, that make it difficult to transfer techniques developed for continuous 
systems to discrete systems, as we will sec later. The construction of the reduced 
space in [s] relics on an isomorphism defined using a connection on the principal 
bundle Q — > Q/G, which is seen as a G-invariant splitting of TQ, that descends to a 
splitting of TQ/G. In the discrete case, we follow the same path, but using what we 
call an affine discrete connection, a minor generalization of the discrete connections 
introduced by M. Leok et al. in [19] which allows us to split {Q x Q)/G. Using, in 
addition, a connection, we are able to derive a reduced variational principle as well 
as reduced equations of motion that split in horizontal and vertical parts. 
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The vertical part of the variational principle and equations turn out to be equiv- 
alent to what is known as the discrete nonholonomic momentum evolution equa- 
tion [10] that determines how the discrete nonholonomic momentum mapping, de- 
fined in terms of the symmetry directions that are compatible with the constraints, 
evolves for discrete mechanical system. The so called horizontal symmetries have 
the property of preserving the discrete nonholonomic momentum. 

Wc specialize the general reduction process in two different settings: the Chap- 
lygin and horizontal systems, where we can go further to obtain reduced systems 
that can be described easily. The first case is already present in the literature and 
we derive the known results from our approach. The other hasn't been considered 
before in the discrete setting to our knowledge; horizontal systems have been con- 
sidered in the continuous setting by Cortes in [11]. Our approach can also be used 
in the case where the configuration space is a Lie Group to re derive the results 
of [23] mentioned above. 

Constraints in a mechanical system on a configuration space Q are usually given 
by two distributions: I) that describes the allowed variations and Ck that limits 
the allowed trajectories. The D'Alembert Principle used to determine the dynamics 
of such systems requires that V = Ck- More general situations, called generalized 
nonholonomic systems, where no connection is made between V and Ck, have been 
considered in [20, 7, 5] . Discrete mechanical systems as introduced in [10] should be 
called generalized in the same sense as above, since the setup for constrained discrete 
mechanical systems consists of V, with the same meaning as in the continuous case, 
and a submanifold Vd C Q x Q which restricts the allowable discrete trajectories, 
and no relation between the two is assumed. 

The layout of the paper is as follows. In Section 2, we sketch very roughly the 
lagrangian reduction of classical mechanical systems, which serves the dual purpose 
of introducing the results in the continuous case as well as the variational approach 
that serves as motivation for our work in the discrete case. Section 3 introduces 
discrete mechanical systems and their symmetries. In Section 4 we introduce some 
tools, including the affine discrete connections, that will be useful in the analysis of 
the reduction process to be carried out in Section 5 and condensed as Theorem 5.11. 
The statement of this result is written in terms that are not obviously defined 
on the reduced system; the purpose of Section 6 is to give an intrinsic version 
of that result, which we achieve with Corollary 6.5. The reconstruction of the 
original dynamics starting from that of the reduced system is explored in Section 7. 
Finally, in Section 8 we study the discrete nonholonomic momentum mapping, 
showing that the vertical part of the reduced variational principle is equivalent to 
the discrete nonholonomic momentum equation. The rest of the paper deals with 
the specialization of the general theory developed so far to particular situations. 
In Section 9 we obtain the equations of motion for systems where Q Q/G is a 
trivial principal bundle. In Sections 10 and 11 we study systems with Chaplygin 
and horizontal symmetries respectively; in the Chaplygin case we obtain intrinsic 
versions that specialize to the results of [11], while in the horizontal case we find 
the discrete version of a type of symmetry whose continuous counterpart had been 
studied in [11]. 

Last, we wish to thank Hernan Cendra for his interest and valuable comments 
on this work. 
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2. Reduction of classical mechanical systems 



In this section wc recall some basic facts of the lagrangian reduction theory of 
(generalized) nonholoiiomic mechanical systems in the presence of symmetry. We 
refer to [3] and [()] for further details. 

2.1. Generalized nonholonomic mechanical systems. 

Definition 2.1. A generalized nonholonomic mechanical system is a quadruple 
{Q, L, TD, Ck) where Q is a differentiable manifold, the configuration space, L : 
TQ — R is a smooth function on the tangent bundle of Q, the lagrangian, T) 
is a subbundle of TQ, the variational constraints or virtual displacements, and 
Ck C TQ is a submanifold, the kinematic constraints. 

For every such system, the action functional is defined by 



where q : [to, ti] — > Q is a smooth curve in Q and q : [to, ti] ^ TQ is its velocity. An 
infinitesimal variation of q is a smooth curve 5q : [to,ti] — >■ TQ . An infinitesimal 
variation is said to have vanishing end points if Sq{to) = and 6q(ti) ~ 0. 

The dynamics of a generalized nonholonomic mechanical system is determined 
by the following Principle. 

Definition 2.2 (Lagrange-D'Alembert Principle). A trajectory of (Q, L,'D,Ck) 
is a curve q : [to, ti] — Q which 

• satisfies the kinematic constraints: {q{t),q{t)) G Ck for all t g [to,ti] and 

• is a critical point of S for the admissible variations: dS(q)(6q) = for all 
infinitesimal variations Sq of q with vanishing end points and such that 
Sq{t) e for aU t e [to, ^i]- 

As usual, the Lagrange-D'Alembert Principle gives rise to a set of equations 
called the generalized Lagrange-D'Alembert equations for {Q,L,T>,Ck)- These 
equations of motion can be written in a coordinate-free way if an affine connection 
V on Q is chosen. Assuming this additional datum, [6] proves the following result. 

Theorem 2.3. A smooth curve q in Q is a trajectory of (Q, L, T>, Ck) if cind only 
if iqit),q{t)) € Ck for all t € [to,ti] and 



(1) - — FL(<z(t),g(t))+BL(g(t),q(t)) G (P,(t))° for all te[to,ti]. 



where ¥L and ML denote the fiber and base derivatives of L, as defined in [6, 
Definition 3] and (23q(t))° C 2^g*(t)'9 the annihilator ofT>g(^fj C T^^i^Q. 

Condition (1) is known as the (generalized) Lagange-D'Alembert equation. 

2.2. Symmetric generalized nonholonomic mechanical systems. In what 
follows, we assume that G is a Lie group that acts (on the left) on Q freely and 
properly, so that the quotient tt : Q Q/G is a. principal bundle with structure 
group C. This action will be denoted by lf{q) for all 5 G G and q € Q; there is 
an induced action of G on TQ given by lj'^{vq) := dlf{q){vq) for a\\ g G G and 
Vq G TqQ, that is called the lifted action. 




to 
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Definition 2.4. G is a symmetry group of (Q, L, P, Ck) if, in addition to tlie 
general assumptions stated above, L, T> and Ck arc G-invariant, that is, if Lo/JQ = 
L, lJ^{V) C 2? and IJ'^{Ck) C Ga' for aU g e G. 

The vertical bundle over Q is a subbundle of TQ with fibers := Tg{l^{{q})). 
We assume that the distribution S with fibers Sq :~ n 2?^ has locally constant 
rank, so that it is a subbundle of TQ. This condition, which applies to many inter- 
esting examples, is verified when, for instance, TqQ = Vq +Vq for all 5 G Q — the 
"dimension assumption" — holds, as is usually considered in much of the literature. 

Assume that G-invariant subbundles of TQ, Ti, U and W can be chosen in such 
a way that Ti. and U are direct complements of 5 in P and V*^ respectively, and W 
is a direct complement of 2? + V*^ in TQ. All together, we have the decomposition 

(2) TQ^W®U®S®U, 

where V — S ® % and V*^ = S ®U. One way to construct the complementary 
bundles "H, U and W is as orthogonal complements for some G-invariant inner 
product on TQ. This type of inner product is usually available as the kinetic 
energy of symmetric mechanical systems. 

Following [S] and [(i], we associate a connection to the decomposition (2). 

Definition 2.5. The unique connection A on the principal bundle tt : Q ^ Q/G 
whose horizontal space is Hor^ = W © H is called the generalized nonholonomic 
connection associated to the system {Q, L,V,Ck) and the splitting (2). 

2.3. Reduction of symmetry. When G is a symmetry group of (Q, L, T>, Ck) we 

define the reduced lagrangian £ : TQ /G — M by t{[{q^ 'z)]g) '■= L{q, q). 

In order to establish a reduced variational principle and study the reduced equa- 
tions of motion [9] and [S] find it more convenient to work on a difFeomorphic model 
of TQ/G as follows. Given a principal connection A an. n : Q Q/G there is an 
isomorphism of vector bundles over Q/G 

aA-.TQ/G ^T{Q/G)(BB, 

where g is the Lie algebra of G, g is the adjoint vector bundle (see [17]) of tt and 
C(a{[<},q]g) ■= dTT{q,q) © [{q, A{q,q))]G- The connection of choice used to define 
aj[ is the generalized nonholonomic connection A introduced above. The reduced 
lagrangian £ can be transported to this new model space as L : T{Q/G) © g -^^ R, 
defined by 

L{x,x,v) := £({aAy^(x,x,v)) = L{q,q), 

where = dn{q,q) and v — [{q, A{q,q))]G- 

The reduced dynamics of the system is defined in [li] by the following principle. 

Definition 2.6 (Generalized Lagrange-D'Alembert-Poincare Principle). A trajec- 
tory of the reduced system determined by the symmetry group G of (Q, L, I?, Ck) 
is a curve ^ := (x, x) (Bv : [to, ^i] — >■ T{Q/G) © g which 

• satisfies the reduced kinematic constraints: fi(t) g Ck ■— Q:a{Ck /G) for 
all t G [to,ti]. 

• is a critical point of the reduced action 




L{x{t),x{t),v{t))dt 
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for some infinitesimal variations Sfi = bx^^-^v of such that, if q is a hft of 
X to Q, &x{t) e d7r(g(t))(X>,(t)) and 8^v{t) e := a^(5/G) with 

aU variations vanishing at the end points. The precise description of the set 
of infinitesimal variations where the criticality condition applies requires a 
special form for the b^v involving the curvature of A (we refer to [6] for 
further details). 

The dynamics generated by this principle can be described using equations of 
motion, called the reduced generalized Lagrange-D'Alembert-Poincare equations, 
which are equivalent to the generalized nonholonomic Lagrange-D'Alembert equa- 
tions of the system (Q, L, D, Ck\ as shown by Theorem 9 of [(>]. 

Theorem 2.7. Let G be a symmetry group of {Q, L,'D,Ck), A be the generalized 
nonholonomic connection on n : Q ^ Q/G associated to a splitting of TQ and 
q : [to,ti] Q a curve on Q. Then the following statements are equivalent. 

• The curve q satisfies {q{t),q(t)) S Gk for all t G [tQ,ti] and (1) holds. 

• The curve /i : [to,ti] T{Q/G)®q given by fj, := {x,x)(Bv — a_A{[{q, q)]G) 
satisfies nit) € Gk and 

n r)f r)f 

for all t G [toi^i] , where B is the reduced curvature of A. 

3. Discrete mechanical systems and symmetries 

In this section we review the notion of discrete mechanical system with nonholo- 
nomic constraints, that is the discrete time analogue of the (generalized) nonholo- 
nomic mechanical systems considered in Section 2. We also consider symmetries of 
such systems. 

3.1. Discrete mechanical systems. 

Definition 3.1. A nonholonomic discrete mechanical system consists of a quadru- 
ple {Q, Ld,T),T>d) where Q and T> are as in Definition 2.1, : Q x Q ^ K. is a 
smooth map, the discrete lagrangian, and 2?^ C Q x Q is a submanifold, the discrete 
kinematic constraints. 

Remark 3.2. The discrete mechanical systems defined above are slightly more 
general than those considered in [10] because we are not requiring that the diagonal 
of Q X Q be contained in T>d. Still, in order to construct a dynamical system, we 
will assume that 2?^ contains the graph of a smooth map Q ^ Q, with the case 
of [1(J] corresponding to the identity. 

Following [10], discrete mechanical systems define discrete dynamical systems 
using a discrete Lagrange-D'Alembert Principle, roughly saying that trajectories 
of the dynamical system are critical points of the discrete action functional 

N-l 

Sd{q-) := ^ Ld{qk,qk+i) 

k=0 

that satisfy the constraints. The following definition makes this notion precise. 
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Definition 3.3 (Discrete Lagrange-D'Alembert Principle). A discrete curve in Q 
is a map q. : {0, 1, . . . , N} — > Q and a variation of a discrete curve q. consists of a 
map (Jg. : {0, 1, . . . N} — TQ such that Sqk S TgfcQ for aU k. A variation is said to 
have vanishing end points if Sqo = and 6q]\r = 0. A trajectory of {Q, Ld, 2?, is 
a discrete curve g. which 

• satisfies the kinematic constraints {qk,qk+i) G T^d for all k and 

• is a critical point of S'd for all admissible variations Jg. ofq.: dSd{q ){Sq.) ~ 
for all infinitesimal variations with vanishing end points Sqk € "Dq^ for all 
k. 

The discrete Lagrange-D'Alembert Principle leads to a set of equations. Indeed, 
if Dj denotes differentiation with respect to the j-th component of a Cartesian 
product, 

EAf-l 
dLd{qk,qk+l){Sqk,Sqk+i) 
k—O 

= Z^fc=i {DiLd{qk,qk+i) + D2Ld{qk-i,qk)){5qk) 
+ DiLd{qa,qi){5qQ) + D2Ld{qN-i,qN){5qN), 
so that q. is a trajectory of {Q, Ld, T), Vd) if and only if, for all fc, 
(3) DiLd{qk,qk-i) + D2Ld{qk-i,qk) e and {qk,qk+i) e 

or, if Xniqk, qk+i) ~ are local equations defining 2?^, 

EM 
Aa,fc w''((3'fe) and Xii{qk,qk+i) = 0, 
a—l 

for some constants Xa,k G and 2?°^ = {^^{qk), . . . , aj*^(gfc)). 

It can be shown that under sufficient regularity of the data, equation (3) has 
solutions, that are unique given sufficiently closely spaced initial data (see [23, 
Proposition 3] and [10]). 

Sometimes it is necessary to consider systems that are forced. Such will be the 
case below when we consider reduced systems, even when the unreduced one is not 
forced. Thus, it is convenient to add the following notion. 

Definition 3.4. A forced discrete mechanical system consists of a discrete me- 
chanical system (Q, Ld, T>, Vd) together with a 1-form fd on Q x Q. We often write 
fd{qo,qi){Sqo,Sqi) ■= fdi(lo,(li){^Qo) + fdilo>qi){^Qi) were /j" : pI{TQ) R and 
fj^ : p\{TQ) -> R, where pj : Q x Q ^ Q denotes the projection on the j-th 
component and p*{TQ) is the pullback vector bundle (see [J]) over Q x Q . 

The dynamics of a forced discrete mechanical system is given by the appropri- 
ately modified discrete Lagrange-D'Alembert principle as follows. 

Definition 3.5. A discrete curve q. is a trajectory of the forced discrete lagrangian 
mechanical system {Q,Ld,'D,'Dd, fd) if 

• the curve satisfies the kinematic constraints: {qk,qk+i) G 'Dd for all k and 

• for all variations 5q. of q. that have vanishing end points and Sqk G Vq^, , 

N-l 

dSd{q.){Sq.) + ^ fd{qk, qk+i){Sqk,Sqk+i) = 0. 

fc=0 
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Just as we found a characterization in terms of equations in the unforced situa- 
tion, q. is a trajectory of a forced system if and only if, for aU k, 

DiLdiqk,qk+i) + D2Ld{qk-i,qk) + fd i'lk,qk+i) + fd{qk-i,qk) e 1^°^ 

(qk.qk+i) e Vd- 

3.2. Symmetric discrete mechanical systems. The main subject of this paper 
is the analysis of discrete mechanical systems with continuous symmetry groups. 
Here we continue under the assumption that the Lie group G acts on Q by making 
the quotient map tt : Q — > Q/G a principal bundle. In addition, we consider the 
diagonal action of G on Q x Q defined by if^'^iqojqi) '■= {l^ (qo) , {qi)) . 

Definition 3.6. G is a symmetry group of {Q, Ld,'D,T)d) if, in addition to what 
was stated in the previous paragraph, Ld and Dd are invariant by l^^^ ^ and P is 
G- invariant by f"'^ . If the system is forced, we also require that the discrete force 
fd be G-equivariant. 

Several structures introduced in Section 2.2 remain useful in the current context. 
In particular, we have the subbundles 5, V*^, H,U,W C TQ and the decomposition 
of TQ given by (2). We also have a nonholonomic connection A associated to that 
decomposition. We denote the horizontal lift associated with A hy h : Q ^q/g 
T(Q/G) — > TQ, that is h'^{'WTr{q)) ~ Vq if G Horj({q) and d7r{q){vq) ~ w-n{q)- The 
following result is straightforward. 

Lemma 3.7. Considering the G-actions given by the lifted action and hy 
^Qxq/gT'IQ/G)^^^ ^^^^^^ (^^(9): iir(g))j ihe map h is G-equivariant. 

For each q ^ Q, the elements ^ G such that ^q((?) G T>q form a subspace in g, 
which may depend on q. For that reason, it is convenient to define the space 

■■^{iq,OeQxg:^Qiq)eVq}. 

Using the projection on the first variable, g^ is a vector bundle on Q. It is easy to 
sec that the differential in the group direction of l'^ establishes an isomorphism of 
vector bundles ^ S. 

4. Some discrete tools 

An important step in the reduction process for classical systems is the passage 
from TQ/G to a model space T{Q/G)(B3, as seen in Section 2. This is achieved using 
the nonholonomic connection. In the discrete case, we follow the same philosophy 
but connections are not the right tool for the task. In this Section we introduce the 
notion of affine discrete connection and use it to construct an isomorphism that, 
eventually, will play the role that a a played in Section 2.3. 

4.1. Afflne discrete connections. 

Definition 4.1. The discrete vertical bundle for the l'^ action of G is the subman- 
ifold 

:= {(qj^iq)) E Q x Q : q e Q, g E G}. 
For g G Q we define V^{q) := H {{q} xQ)cQxQ. 
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Definition 4.2. The composition of vertical and arbitrary elements of Q x Q with 
the same first element is defined by • : xq{QxQ)^QxQ with 

{Qo,lf{qo)) ■ {qo,qi) ■= {qoJgiqi))- 

We denote by Z*^ the conjugation action of G on itself, that is, lf{h) :— ghg^^ 
for all g,h ^ G. 

Definition 4.3. Let 7 : Q — > G be a smooth G-equivariant map with respect 
to l'^ and l'^, r {(g, '^(^^((z)) : q € Q} and Hor C Q x Q be a G-invariant 
submanifold such that F C Hor. For each q € Q let Hor{q) Hor n {{q} x Q) 
and Hor^{q) := p2{Hor{q)). 

We say that Hor defines the affine discrete connection Ad on the principal bundle 
TT : Q Q/G if, for each q e Q and aU qi e lci{q}) n Hor'^{q), 

(4) Hor^{q) C Q is a submanifold and Tq,Q = Tq, {l^Hq})) © Tq,Hor^{q) 
(see Figure 1(a)). We denote Hor by HorA^, and call 7 (or even F) the level of 

Remark 4.4. Recalling the notion of transversality of submanifolds (see [14]), 
condition (4) is equivalent to requiring that for each q £ Q, 

dim{Hor^{q))^dimQ-dimG and l2i{q})(nHor^{q). 

Remark 4.5. The notion of discrete connection introduced in [ I !)], coincides with 
that of an affine discrete connection where the level F is the diagonal of Q x Q. 
Since the diagonal oiQxQ plays the role of a "null element" for composition, affine 
discrete connections need not contain the "null element" in their horizontal space, 
just like affine spaces need not contain the null element of a vector space. 

Remark 4.6. The idea behind the introduction of the previous definitions is that 
QxQ should be a discrete version of TQ. Even though this is a powerful idea, there 
are several important differences between those spaces. Using the projection on the 
first factor and the standard projection, both spaces are fibcred over Q, but TQ is 
a vector bundle, while Q x Q is usually not one. In particular, tangent vectors at 
the same base point can be added, whereas there is nothing similar for elements 
of Q X Q. A partial fix for this problem is the composition operation introduced 
in Definition 4.2: it provides a way of combining vertical elements of Q x Q with 
arbitrary elements of Q x Q based at the same point (that is, with the same first 
component). Even though this is not a complete analogue of addition, it is enough 
to handle discrete connections. 

Connections are an important tool in differential geometry. Essentially, for prin- 
cipal bundles n : Q ^ Q/G, they provide compatible splittings 

(5) TgQ = ® Hor{q) for aU q e Q. 

In the discrete setting, we want to be able to split {q} x Q in vertical and (some) 
complementary space. Not having an addition operation we choose a more geomet- 
ric view of the problem and define affine discrete connections in terms of spaces 
that are complementary to in the precise sense of Definition 4.3. Furthermore, 
this definition is equivalent to being able to decompose {q} x Q = V^(g) • Hor{q) 
(at least in a neighborhood of F), which is the composition-analogue of (5). 



10 



JAVIER FERNANDEZ, CORA TORI AND MARCELA ZUCCALLI 



Proposition 4.7. Let Ad be an affine discrete connection of level 7. Then, there 
exists U , a G-invariant open neighborhood ofT in QxQ, such that, for all (qo, qi) S 
Q X Q with 7r(gi) sufficiently close to 7r(qo)j there is a unique g G G such that 

(6) (go,9i) = (go,Z?(9o))-(go,/^-i(gi)), 

with (goJ^-i(gi)) e HorAAlo) n U. 

Proof Since ?^(,„)(go) e /g({go}) niJor^^ (go), the condition ;g({go}) HT Hor^^ (go) 
implies that Hor\ (go) intersects all orbits of G that are close to the one through 
go (see Figure 1(b)). Therefore, for any gi with 7r(gi) sufficiently close to 7r(go), 
there are g £ G such that /^_i(gi) G Hor^^qo). Furthermore, again by the 
transversality condition, there is an open set Vq^ C Q where the intersection 
Z^({gi}) n Hor^^qo) fl Vq^ consists of a single point. Hence, there is a unique 
g £ G such that ^^_i(gi) G i/or^^ (go) H V^,,. Finally, taking into account the 
smoothness of Hor, U is constructed by gluing the sets Vqg . □ 




Figure 1. G orbit and Hor^^ (go) in Q 



Remark 4.8. The existence of a level function as a datum for an affine discrete 
connection may seem a bit strange at first. However, if a decomposition like (6) is 
expected, by taking gi = go, we see that (go, '^-i (<?o)) G HorA^ (go), so 7(go) := g~^ 
is a level function for Ad- 

Definition 4.9. Given an affine discrete connection Ad we define its discrete con- 
nection 1-form Ad ■ Q X Q ^ G hy Ad{qo,qi) g where g is the element of G 
that appears in decomposition (6). 

Remark 4.10. Notice that we will use the same letter to name both the connection 
Ad and its discrete connection 1-form. Also, calling Ad a 1-form can be misleading 
since it is not a 1-form in the usual sense but, rather, a function. Still, the name is 
coming from the fact that when Leok et al. introduced it in [19] they also introduce 
the notion of "discrete /c-form" , and Ad is a discrete 1-form in that sense. 
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Remark 4.11. Due to Proposition 4.7, Ad is only defined in a G-invariant open 
set U C Q X Q. In what follows, we will abuse the notation and pretend that Ad 
is defined everywhere not to make the notation too cumbersome. 

Proposition 4.12. Let Ad be an ajfine discrete connection on the principal bundle 
TT : Q ^ Q/G. Then, for all {qo, gi) G Q x Q and go,gi G G, 

(7) Adil'^,iqo)J^M)) ^ 9iMQo,qi)go\ 

as long as both sides are defined. Conversely, given a smooth function A : Q x Q ^ 
G such that (7) holds (with Ad replaced by A), then Hor := {((i'o,'Zi) (z Q x Q : 
■^(^Oj^i) = e} defines an affine discrete connection Ad with level j{q) := A{q,q)~^ 
and whose discrete connection 1-form is A. 

Proof Writing Ad{qo,qi) = .9 and Ad{lf„{qo),lg,{qi)) = .9, by definition, 
(90,91) = (90,^^(90))- (90, /^-i (91)) 

il^,{qo)jfA<li)) = il^iqo)jfil'^M))) ■ a?o(9o),^f-ia?,(9i))) . 

" V ' 

eHor^Jlgiqo)) 

By the G-invariance of Horj^^, since (go, ^^-1 (91)) £ Horj^^{qQ), it follows that 
(^^c(9o),^L-i(9i)) e Horj,,{ll{q^)). We can write 

EHorAjlfoilo)) 

Identity (7) follows, then, from the uniqueness of the decomposition (6). 

The converse result is, mostly, routine checking and we will omit the details. 
That Hor is a submanifold follows from the easily proved fact that e G G is a 
regular value of A. The same proof shows that p2{Hor{qo)) is a submanifold of Q 
of the right dimension. The G-equivariance of A shows that Hor is G-invariant. 
An application of (7) leads to establishing the transversality condition for Hor. 
Finally, taking 7(9) :— A{q,q)~^, one concludes that Hor defines an affine discrete 
connection of level 7. □ 

Affine discrete connections have discrete horizontal lifts, just as regular connec- 
tions do. 

Definition 4.13. Let Ad be an affine discrete connection on the principal bundle 
TT : Q — > Q/G. The discrete horizontal lift h^ : Q x Q / G ^ Q x Q is given by 

^d°(^i) (go,9i) ^ (90,91) e -H^or^d and 7r(gi) = ri. 

We deimehf ■.= p2ohf . 

Remark 4.14. The map hf{ri) is well defined, provided that ri is sufficiently 
close to 7r(go), since for any qi G 7r~^(ri), 

(8) '^*(^i)^=(gO''X(,o,,0-(«i))- 

Lemma 4.15. Let Ad be an affine discrete connection on tt : Q Q/G and 

consider the natural actions of G on Q x {Q/G) (extending trivially on the 
second factor) and the diagonal action on Q x Q. Then, h^ is G-equivariant. 
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Proof. It is a direct computation using (8). □ 

For completeness we mention that having a smooth and C?-equivariant map rj : 
Q X {Q/G) Q X Q such that n op2 o i]'^ = id\Q/G for all g G Q is equivalent to 
having an affinc discrete connection whose discrete horizontal lift is rj. 

Remark 4.16. Discrete connection forms Ad and discrete horizontal lifts may not 
be defined everywhere. Indeed, li : Q x [Q/G) — >■ Q x Q is defined everywhere, 
then for any q € Q, the map r i— !■ h1i{r) is a global section of the principal bundle 
TT : Q ^ Q/G, so that the bundle is trivial. Hence, for nontrivial principal bundles 
/i^ and, consequently. Ad can only be defined in some open set of the total space. 

On the other hand, if Q has a G-invariant riemannian metric — as it usually 
happens for mechanical systems — the following construction provides a discrete 
affine connection. Let 7 be a level function and define 

i7or := {(go,cxp^(^^)(wi)) G Q X Q for some vi e (V^^^^^)^ <Z T^(^q,)Q}, 

where exp and _L are those of the riemannian metric and 7(90) ^^(qo) ('^o)' 
can be checked that Hor defines a discrete affine connection of level 7. Notice that 
the domain of the associated discrete affine connection form and discrete horizontal 
lifts are limited by the domain of the exponential mapping. 

Example 4.17. Consider Q := with the action of G := R given by l^{q) :~ 
{x,y + 5), where q = {x,y). Clearly, pi = tt : Q — > Q/G is a (trivial) principal 
bundle with structure group G. In order to define an affine discrete connection on 
this bundle we need to find a G-invariant manifold Hor that is "complementary" to 
= {(go, qi) & Q X Q : xo = xi} in QxQ. We consider complements to l%{{qa\) 
for qo E Q. A G-invariant family of curves that is complementary to the orbits is 
Hor^^\qo) = {qi & Q ■ yi — ya ^ b{xi + xo){xi — xo)/2} for a parameter 6 e M. 
Since ^^({go}) DHor^'^^qo) = {qo}, the only thing we have to prove in order to see 
that Hor = {{qo, qi) ^ Q x Q : yi — yo = b{xi + xo){xi ~ xo)/2} defines an aSine 
discrete connection is that Tq„Q = TyJ^^qo}) © Tq„Hor^'^\qo), which is evident 
because TqJ^iUo}) = {dy\qj and Tq^^H or^'^\qo) = {dx\^^^ + bxody\^^). We denote 
the discrete connection defined by Hor with 
Since 

(90, gi) = {qo,l^{qo)) ■ {qo,l'^-i{qi)) for g = yi-yo- ^6(.ti +xq){xi - xq), 

d 

A/'d{qQ,qi) = yi - yo - b{xi + xa){xi - a;o)/2 
h''d{ri) = {qo, {ri,yo + b{ri + xo){ri - xo)/2)). 

The fact that A''^ and are defined everywhere is compatible with Q ^ Q/G 
being a trivial principal bundle. 

4.2. Isomorphisms associated to an affine discrete connection. When work- 
ing with symmetric discrete mechanical systems, one is led to consider the space 
{Q xQ)/G. It is convenient to have a different model for this space. In this section 
we construct such a model associated to an affine discrete connection Ad- 



we have 
(9) 
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We start with a special case of a general construction called the associated bundle 
of a principal bundle (see Chap. 1, Sect. 5 of [17]). 

Definition 4.18. Let tt : Q — > Q/G be a principal bundle and consider the action of 
G on Q X G defined by if^'^iq, w) := {lf{q), lf{,w)), with lf{w) := gwg~^. Being tt 
a principal bundle, the quotient G :^ {Q x G)/G by this action is a manifold, called 
the conjugate associated bundle. The quotient map is denoted by p : Q x G — G. 
The projections onto each of the two components of Q x G induce smooth maps 
pQ/^ : G ^ Q/G and p^/^ : G ^ G/G. The first, p^/^, turns G into a bundle 
over Q/G with fiber G. 

It is convenient to define 

Fi :gxGx (g/G)^g with Fi{qo,wo,n):=l2„{Win))- 

Extending the G action l'-^^^ to Q x G x [Q /G) trivially on the last component 
and considering the G action on Q, a simple computation using the G-equivariance 
of /i^ shows that Fi is G-equivariant. 

Proposition 4.19. Given an ajfine discrete connection Ad on n : Q Q/G, let 
^Aa -QxQ^QxGx {Q/G) and : Q x G x {Q/G) Q x Q be defined by 

^.Adlgcgi) := (go, A('7o,'?i),7r((7i)) 
^.Ad(9o,Wo,ri) := {qo,Fi{qo,wo,ri)). 

Then, and Ad o,^^ smooth and mutually inverses (when restricted to the 

domain of Ad)- Furthermore, considering the diagonal action of G on Q x Q and 
^QxGx(Q/G)^^^^ ^^-j ._ QQ^q^'^^iG^yj^'j j.^^^ both maps are G-equivariant, so that 

they induce diffeomorphisms ■ {QxQ)/G ^ Gx{Q/G) and Ad ■ Gx{Q/G) 
{Q X Q)/G. 

Proof. That $_4^ and ^ Ad ^-rc smooth is clear from the definition. Checking that 
they are inverses and the G-cquivariance is done by direct computation. □ 

Remark 4.20. In fact, when Ad and are not globally defined the maps 
and Ad are diffeomorphisms between open neighborhoods of of F/G and {Q x 
{e})/G X {Q/G). 

The commutative diagram (10), where tt is the quotient map and T := o tt, 
shows some of the spaces and maps we have introduced and that will be used later, 
when analyzing the behavior of symmetric discrete mechanical systems. 

(10) QxQ -l^QxGx{Q/G) 




pxid 



{Q^Q)/G—^Gx{Q/G) 

Lemma 4.21. Let {qo,qi) <E Q x Q and {Sqo,Sqi) G T(^q„,qi}{Q x Q). Then, 
dT{qo, qi){Sqo, 5qi) = {dp{qo, Ad{qo, qi)){Sqo, dAd{qo, qi){Sqo, Sqi)),dTT{qi){5qi)). 



Proof. Compute dT, using T = (p x id) o ^Ad- 



□ 
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Example 4.22. In the context of Example 4.17, we have that Q ~ {Q/G) x G with 
the G-action on Q correspondmg to left multiplication on G. Then G ^ (Q/G) xG 
with p{{ro,ho),wo) i— > (ro,wo). The map p has a section s(ro,wo) •= ((^o,0),wo) 
(where we are using implicitly the isomorphism just defined). The isomorphisms 
that appear in Proposition 4.19 are 

*^^('7o,'?i) = {{xo,yo),yi - ijo - b{xi + xo){xi -xq)/2,xi) 

*^6((?'o, ha), Wo, ri) = ((ro, ha), {n, ha + b{ri + ro)(ri - ro)/2 + wq)). 

Remark 4.23. The model spaces Q x G x (Q/G) and G x (Q/G) are by no means 
unique. In fact, given an affine discrete connection Ad it is possible to consider the 
isomorphism aA^ '■ [Q Q) I G ^ {Q / G x Q /G) Xqjg G defined by 

a^d(7r('7o,'7i)) = ((7r((7o), 7r((7i)), p(go, -4d((7o, ^i)))- 

This isomorphism and the corresponding model space are very close to the data 
used in the analysis of the continuous case in Section 2.3 and is introduced by [19] 
in their proposed reduction of discrete unconstrained systems. Still, the presence 
of the fibcrcd product in the model space makes it harder to work with compared 
to the simpler Cartesian product that appears in G x (Q/G). 

5. VARIATIONS AND REDUCED VARIATIONS 

In this section we analyze the relationship between the dynamics of a symmet- 
ric discrete mechanical system on Q x Q, which induces a dynamics on (an open 
neighborhood of Q x {e} x (Q/G) in) (Q x Q)/G and the dynamics on (an open 
neighborhood of [Q x {e})/G x {Q/G) in) the isomorphic model G x {Q/G) of 
{Q X Q)/G. In order to describe the dynamics on this last space we start by 
introducing some relevant notions. 

5.1. Reduced lagrangians. Given {Q,Ld,'D,Vd) with a symmetry group G and 
an afhne discrete connection Ad we construct a dynamical system on Q x G x {Q/G) 
using the isomorphism . In particular, the dynamics will be determined using 
the function Ld '■= Ld o Since Ld is G-invariant and, taking into account the 

symmetry properties of the system, Ld is G-invariant, Ld and Ld induce maps on 
the corresponding quotient spaces. We denote the map induced by Ld on Gx {Q/G) 
by Ld- Therefore, 

Ld{qo,wo,ri) = Ld{qo,Fi{qo,wo,ri)) and Ld{p{qo,'Wo),ri) = Ld{qo,wa,ri). 
The following diagram shows all the relevant maps introduced so far. 




{Q^Q)/G^^Gx{Q/G) 

Associated to the reduced discrete lagrangian Ld we define a reduced discrete action 
Sd by Sd{v.,r.) := J2k ^d{vk,rk+i)- 
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Example 5.1. Consider the discrete mechanical system (Q, LdjVjVd), where 
Q := with points q = {x,y), 
Ld{qo,qi) := m((.Ti - a-o)^ + (yi - 2/o)^)/2, 
V, := {xd,\^ + ydy\^ eT,Q:y = xx} ^{d^^- xdy\^) C T,Q, 
T^d := {(go, gi) e Q X Q : 2/1 - 2/0 = (a;i + xo){xi - xo)/2}, 

that originates as a discretization of a free particle in subject to a nonholo- 
nomic constraint. This system can be readily solved using the discrete Lagrange- 
D'Alembert equations (3). Rather than following that path we will make a detailed 
study of its reduction in the following sections in order to illustrate some of the 
techniques developed so far. 

The group G := M acts on Q as considered in Example 4.17, turning G into a 
symmetry group of {Q,Ld,'D,Vd). Furthermore, for arbitrary 6 G K, defines 
isomorphisms $_^b and <f>_^b that can be used to study the system. Since, by 
Example 4.22, G x (Q/G) ~ {Q/G) x G x (Q/G), the relevant lagrangians are 
Ld{{ro, ho),wo,7'i) = m{{ri - ro)^ + (wq + b{ri + ro)(ri - ro)/2)^)/2 
Ld{ro,WQ,ri) = Lrf((ro, 0), Wq, ''i)- 

Back in the general setting, we notice that, since Ld is G-invariant and Q x 
G X (Q/G) is a Cartesian product (so that differentials on that space decompose 
in terms of differentials on each of the spaces), 

dldivo, ri)((Suo, Sri) =dXd(p(go, wo), ri)(dp(go, wo)iSqo,SwQ),Sri) 

^^^^ ^dLd{qo,'Wo,ri)(Sqo,Swo,Sri) 

=DiLd{qo,wo,ri){Sqo) + D2Ldiqo,wo,ri){Swo) 

+ D3Ld{qo,Wo,ri){Sri). 

In Section 5.2, we study more carefully the second term in this last equality. 

Lemma 5.2. Considering the G actions iQxGx(q/g) ^ iQ ^ ^^^^ their lifted actions 
to the corresponding tangent spaces, the following morphisms of vector bundles over 
QxGx {Q/G) 

DiLd : p*iTQ R and D2Ld : P2TG ^ M. 

are G-equivariant, where pi and p2 are the projections of Q x G x {Q/G) onto the 
corresponding factors. 

Proof. The G-invariance of Ld leads to the G-equivariance of dLd : T{Q x G x 
(Q/G)) M. The statement follows by decomposing T{Q xGx {Q/G)) = plTQ® 
p*2TG®plT{Q/G). □ 

5.2. Mixed curvature and reduced forces. Given a principal bundle tt : Q ^ 
Q/G., a connection A and an affinc discrete connection Ad we construct an ob- 
ject that, somehow, compares the notions of continuous and discrete horizontality 
introduced by A and Ad- 

Definition 5.3. The mixed curvature Bm of A and Ad is the morphism of vector 
bundles over QxQ.B^: T{Q x Q) ^ Ad*{TG) defined by 

Bm{qo,qi){Sqo,6qi) := dAd{qo,qi){HorA{Sqo), HorA{Sqi)), 
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where Horj\[5q) denotes the ^-horizontal part of 5q. Taking into account the 
isomorphism T[Q x Q) ^ Pi{TQ) (Bp2{TQ), the mixed curvature decomposes as 
B.,n = B+ + B;-^, where B+ : p^TQ ^ Ad*{TG) and : plTQ ^ Ad*{TG) arc 

I3m{(la,(li){^n) ■=l3m{qo,qi){0,Sqi) = D2Ad{qo,qi){HorA{6qi)) 
Bmiqa.qi)iSqQ) ■.=B„i{qQ,qi){Sqo,0) = DiAd{qa, qi)iHorAiSqa)). 

Remark 5.4. An obvious difference between the mixed curvature and the curva- 
ture of a connection is that the former can be represented as a 1-form, while the 
latter is represented as a 2-form. 

Remark 5.5. The mixed curvature measures how ^-horizontal deformations at 
(go, qi) depart from the level manifolds of Ad- 

Example 5.6. Wc continue in the context of Example 5.1. In order to compute 
the mixed curvature (and continue analyzing the reduction of the system) wc in- 
troduce a connection on the principal bundle tt : Q — > Q/G. For that purpose, the 
splitting (2) is 




The connection A with horizontal space Horj\, = W © "H has connection 1-form and 
horizontal lift given by 

A{xdx\^^ ydy\^) = y - XX and h'i{dr\^) = dx\^ + rdy\^. 
We can now compute the mixed curvature Bm associated to A and ^J^: 

Bm{qo, qi){SqQ, Sqi) ^dAd{qo, qi){HorA{Sqo), HoTAiSqi)) 

=dw„\ Ab, ^®(dyi\ ~ dyo\ - bxidxi\ + bxodxo\ ) 

i^oidoc,]^^ + xody,\^J + ci{d^,\^^ +xidy,\^^)) 
=(1 -b){cixi -coXo)9u,oL6(,Q,,^)- 

Notice that depending on the "relative slopes" 1 and 6 of ^ and A'^, the mixed 
curvature vanishes or not. 

Just as regular curvature, the mixed curvature can be seen as a "reduced object" 
defined on (Q x Q)/G or, using the isomorphism ^Aa^ on G x (Q/G). We study 
this object next. 

Since Ad : Q x Q ^ G is G-cquivariant for iQ^Q and l"-^, considering the lifted 
actions on T{Q x Q) and TG, dAd induces a morphism dAd on the quotient vector 
bundles: 

T{Q X Q)/G^^TG/G 

(Q X Q)/G—^G/G 

Let a : Gx [Q/G) [Q/G) x (Q/G) be cr(wo, n) := {p'^/'^{vo), rj; then, we have 
the vector bundles a*T{Q/G x Q/G) over G x (Q/G) and {p x id)*a*T{Q/G xQ/G) 
over Q xGx {Q/G). We define x : (p x id)*a*T{Q/G x Q/G) T{Q x Q) by 
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x{qQ,wo,ri){6rQ,Sri) := {h'^" {Sro), h"^^ {Sri)) for qi -.^ Fi{qo,wo,ri). As Fi and the 
horizontal hft are G-equi variant, x descends to a morphism of vector bundles x- 

a*T{Q/G X Q/G) T{Q x Q)/G 



G X (Q/G) ^ (g X Q)/G 

It will be convenient to consider the composition dAd o as a morphism of the 
vector bundles a*T{Q/G x Q/G) {Ad o vE'^^)*(rG/G) over the same base space 
G X {Q/G). Even better, since Ad o ~ pG/G have the following notion. 

Definition 5.7. The reduced mixed curvature Bm is the morphism of bundles over 
G X (Q/G) Brrr-^dAdOx: a*T{Q/G x Q/G) ^ (p^/G opi)*(TG/G). Explicitly, 

^™ («o ,n){Sro, 5n ):^Bra (go , Qi ) {h"" {Sro), h'^' {Sn ) )] 
= [dAd{qo,qi){h''"{Sro),h'J'{Sri))], 

where vq = p{qo,wo), qi := Fi{qQ,wo,ri) and [•] denotes the equivalence class in 
TG/G. Associated to the decomposition B„i ~ i3+ + i3~ there is a decomposition 
B,n ^B++B- where B+ : p^r(Q/G) ^ (p^/^ o pi)*(rG/G) and B- : (p^/^ o 
Pi)*r(Q/G) ^ (p"=^/^ opi)* (TG/G) are defined by 

B+{vo,ri){Sri) ■.^[B+{qo, qi){h'^Hdri))] = [D2Ad{qo, qi){h''HSri))] 
B~{vo,ri){Sro) :=[S„(go, 9i)(/i* (^o))] - [DiAd{qo, qi){h''° {Sro))]- 

Now we return to the discrete lagrangians introduced in section 5.1. In particular, 
we associate a new notion to the second term in the last equality of (11). 

By Lemma 5.2, D2Ld '■ P2TG — > R is a G-equivariant morphism of vector bundles 
over Q X G X {Q/G). Thus, it induces a morphism of the corresponding quotient 

vector bundles D2Ld'. 

{p;tg)/g^^r 

G X {Q/G) 

Notice that {plTG)/G ~ {p^'^ opi)*TG/G. Then, Dz^d o B„, is a well defined 
morphism of vector bundles. 

Definition 5.8. The reduced discrete force is the morphism of vector bundles over 

G X (Q/G) defined by Fd := D2Ld o Bm : cr*(r(Q/G x Q/G)) ^ M (essentially a 
1-form over G x {Q/G)). Concretely, 

(12) Fd{vo,ri){Sro,6ri) i?2^d(go, wo,ri)dAd{qo, gi)(/i* (<5ro), h"^' {Sr^)) 
where vq = p{qo, wq) and qi :— Fi{qo, wq, ri). Once again, using T{Q/G x Q/G) ~ 
pIT{Q/G) (Bp*2T{Q/G), we define F+ Ds^d o B+ and F^ Da^d o i5„, so 
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that Fd = + F^. Explicitly, 

F+{vo,n){Sn) := D2Ld{qo,wo,n)D2Ad{qo,qi){h''H^n)) 
Fd ivo,ri){Sro) := D2Ld{qo,Wo,ri)DiAd{qo,qi){h''° (Sro)). 

Example 5.9. Continuing the analysis of the discrete mechanical system intro- 
duced in Example 5.1, we compute its reduced discrete force with respect to A and 
If {vo,ri) ^ {ro,wo,ri) G [Q/G) x G x (Q/G) we can take q^ := (ro,0) and 
qi := lu<o{h^d {fi)) ~ {ri,Wo + b{r'f — rf^)/2) in the computation. Since 

D2Ld{{rQ,Q),wo,ri) = m(wQ + b{rf - rl)/2)dwo\^^, 

using the results of Example 5.6, we have 

Fdiro, Wo, ri){c„dr„\^^, cidr^l^J = (m{wa + b{rl - ro)/2)dwo|,^J 

((l-6)(ciri ~ coro)dyj„\^J 
= m(l - b){ciri - Coro){wo + b{rl - r^)/2). 

Equivalently, 

Fdiro,wo,ri) = mil - b)(wQ + b{rl - rl)/2){ridri - rodro). 

5.3. Reduced dynamics. In this section we relate the dynamics of a symmetric 
discrete mechanical system to the dynamics of a "reduced" system, defined using a 
variational principle. 

Remark 5.10. The analysis of the symmetric system leads us to consider some 
objects defined on the symmetry group G, that is a Lie group. We review very 
briefly a notation that is convenient and used in this subject. When wqjWi € G 
and Swo G Tw^G we define 

wiSwq dL^„i(wo)((5ti'o) e TujywoG and SwoWi dR.,^^{wo){Swo) S T^„wiG, 

where L^-^ and R^i denote the multiplication by wi on the left and right, respec- 
tively. Analogously, when ao 6 T*^^^G, 

wiao '-^{dL^^-iiiviWa))* {ao) e T*^^,^^G and 
aowi :^{dR^-i{woWi))*{an) G T^^wfi. 

Last, we notice that if ao £ T^^G and (5t«o G T^oG, the following identities hold 

aoiSwo) = {wiao){wi5wo) = {aoWi){Swowi). 

Theorem 5.11. Let q. be a discrete curve in Q, rk TT{qk), Wk Ad{qk,qk+i) 
and Vk := p{qk,Wk) be the corresponding discrete curves in Q/G, G and G. Then, 
given a discrete mechanical system {Q,Ld,'D,'Dd) with symmetry group G, the fol- 
lowing statements are equivalent. 

(1) {(ik,qk+i) £ T^d for all k andq. satisfies the variational principle dSd{q ){Sq.) = 
for all vanishing end points variations dq. such that Sqk S "Dq^, for all k. 

(2) q. satisfies the Lagrange-D'Alembert equations (3) for all k. 
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(3) {vk,rk+i) e Vd := $^,(Pd/G) for all k and dSd{r.,v.){Sr.,6v.) = for 
all variations {dv.,Sr.) with vanishing end points such that 6rk G 1?,.^ 
dTr{qk){'DqJ and 

^^^^ Svk := dp{qk,Wk){h'^''{Srk),dAd{qk,qk+i){h'"'{5rk),h'"'^^{Srk+i))) 

+ dp{qk,Wk){{(,k)Qiqk), dAdiqk, qk+i){{^k)Q{qk), {£.k+i)Q{qq+i))), 

where (qkAk) G 9^ ■ 

(4) {vk,rk+i) € T^d for all k and {v.,r.) satisfies the following conditions for 
each fixed {vk-i,rk,Vk,rk+i). 

• e T* (Q/G) defined by 



(15) 



(j) ■.=DiLdiqk,Wk,rk+i) o h'"' + D3Ld{qk-i,Wk-i,rk) 
+ F^{vk,rk+i) + F+{vk-i,rk) 

vanishes on 'Dj.^, i.e., 

(16) </'eP°^. 

• Tp ^ Q* defined by 

(17) ■■=D2Ld[qk-i,Wk-i,rk)w]:l^ - D2Ld{qk,Wk,rk+i)w^^ 



vanishes on g^, i.e 



V 

, I.e., 

V \o 



(18) V'e(0,"J 

Proof. First wc study the equivalence of conditions (g/c, Qfe+i) G I?d and (ufc, r/j+i) G 
T>d- If iqk,qk+i) G Vd, it follows that Tr{qk,qk+i) G Vd/G and (i;j,,rfe+i) = 
*^d(^(9fc,9fc+i)) e ^AaC^d/G). Conversely, if{vk,rk+i) = ^Aai^iqk, qk+ij) G Vd, 
we have that Tf{qk,qk+i) G Vd/G, so that {qk,qk+i) G Pd- 

Next we tackle the equivalence between the variational principles and equations. 

1 <^=> 2. Is a standard computation using calculus of variations. 

3^1. Given a variation Sq. as in the statement, we let Srk '■= d-K{qk){5qk) and 
6vk using (14) for ^fe := A{qk){6qk), so that 5qk = h'^''{Srk) + (6)Q('7fc)- Using 
Lemma 4.21 to compute dT, 

dSd{q.){Sq.) =2^^^^ dLd{qk,qk+i){Sqk,Sqk+i) 

(19) =V^ ^ dL,(T(gfc,g,+i))dT(%,Jg,+i) 

,JV-1 



dL<j(wfe,rfe+i)((5w/c,(5?'/c+i) = fiS'rf(r., «.)(()>., (5i;.) = 0, 



'fe=0 

where the last equality holds because (r. , u.) satisfies condition 3. Hence, q. satisfies 
condition 1. 

1 => 3. Given a variation ((5r.,(5u.) as in the statement, let 8qk := h'^''(Srk) + 
i^k)Qiqk)- Then, using Lemma 4.21 and the explicit form of (14), wc have that 
{Svk, 6rk+i) = dT{qk, qk+i){Sqk, Sqk+i)- A computation similar to (19) shows that 

dSd{v.,r.){5v.,Sr.) ^ dSd{q.)i6q-) ^ 0- 
since the variation dq. satisfies condition 1. Hence, condition 3 holds. 
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3 4. For variations {Sv.,5r.), writing Svk = dp{qk,Wk){6qk,Swk), 



dSd{v.,r.){Sv.,Sr.) = 2^ dLd{vk,rk+i){6vk,Srk+i) 



N-l 

N-l 
k=0 



dLd{qk,Wk,rk+i){5qk,5wk,5rk+i) 
(DiLd{qk,Wk,rk+i){dqk) + D2Ld{qk,Wk,rk+i){dwk) 



+ D3Ld{qk,Wk,rk+i){Srk+i))- 

Using A to decompose Sqk = h'^''{Srk) + {£,k)Q{qk) and rearranging indexes, for 
vanishing end point variations, we obtain 



(20) 



N-l 



dSd{v.,r.){Sv.,Sr.) = ^ (DiLd{qk,Wk,rk+i) o h'^'' + D3Ld{qk~i,Wk-i,rk) 



k=0 



+ D2Ldiqk,Wk,rk+i)DiAdiqk,qk+i) ° h'"' 
+ D2Ld{qk-i,Wk-i,rk)D2Ad{qk-i,qk) o h'^'' )(frfc) 



^ {DiLd{qk,Wk,rk+i)({^k)Q{qk)) 



k=0 

+ D2Ld{qk,Wk,rk+i)dAd{qk,qk+i){{£,k)Q{qk), {^k+i)Q{qk+i))) ■ 

Since the variations 6r. are independent from those generated by the ^. , we conclude 
that condition dSd{r.,v.){Sr., 6v.) = from point 3 in the statement is equivalent 
to the vanishing of the first and second summations in (20) independently, for all 
vanishing end points variations 6r. with Svk S I^r^ for all k and for all ^. with 
(9fe, Cfe) e for all k and = d = 0. 

Recalling that the variations 6r. are independent and (13), it is clear that the 
vanishing of the first summation in (20) is equivalent to the fact that 0, defined 
in (15); satisfies condition (16). 

Before we prove the equivalence of condition (18) and the vanishing of the second 
summation in (20) we need two auxiliary computations. On the one hand. 



dAdiqk,qk+i){{^k)Qiqk),{(k+i)Qiqk+i)) 
- A. 

~ di 



cxp{t£,k+i)Ad{qk, qk+i) exp{-t^k)) 



t=o 



= ^k+lWk - Wk^k, 

where the notation used in the last equality is the one introduced in Remark 5.10. 
On the other hand, since Ld{lf{qk),Wk,rk+i) = Ld{qkjg-iiwk),rk+i), 



DiLd(qk,Wk,rk+i){{^k)Q{qk)) 



dt 



t=o 



Ld{l%f^t^^^{qk):Wk,rk+i) 



D2Ld{qk,Wk,rk 



+1. 



dt 



{cxp{-t^k)wk cxp(t^fe)) 



t=o 



= D2Ld{qk,Wk,rk+i){-ikWk +Wk£,k)- 
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Using the previous computations we see that the vanishing of the second summation 
in (20) is equivalent to 

=2^^_-^ [D2Ld{qk,Wk,rk+i){^k+iWk - ^fcWfc)) 

^Z^fc-i {D2Ld{qk-i,Wk-i,rk)wi^_^ ~ D2Ld{qk,Wk,rk+i)w^. )(Cfc) 

for all ^. with {qk,^k) ^ Q^- It follows immediately that this last condition is 
equivalent to the fact that ip, defined in (17), meets condition (18). □ 

From the proof of Theorem 5.11 we isolate the following partial result. 

Lemma 5.12. With the notation of Theorem 5.11, the following assertions are 
equivalent. 

(1) The discrete curve q. satisfies dSd{q ){5q'^) = for all vanishing end points 
variations such that 5qk £ Sq^. for all k. 

(2) The curve {v.,r.) satisfies condition (18) for ip defined by (17). 

Remark 5.13. We notice that the horizontal equations that appear in item 4 
of Theorem 5.11 contain force terms that, because of (12) are a composition of 
derivatives of the reduced lagrangian and a term involving A and Ad- Since in 
the continuous setting, the forces that appear in the analogous equation involve 
the reduced curvature of the nonholonomic connection, we chose to call the corre- 
sponding term in the discrete setting the reduced mixed curvature. In this case, we 
recover the standard result that the vanishing of the mixed curvature implies the 
vanishing of the discrete forces in the reduced system. Other than that, we do not 
have a good reason to call Bm or Bm "curvatures" . 

Remark 5.14. For classical mechanical systems, the Lagrange-D'Alembert equa- 
tions (1) arc second order differential equations, while their reduced counterparts 
that appear in Theorem 2.7 are second order in the Q/G variables but only first 
order in the g ones. Analogously, for discrete mechanical systems, the discrete 
Lagrange-D'Alembert equations (3) are recurrence equations of second order while 
the reduced equations obtained in item 4 of Theorem 5.11 are second order in the 
Q/G variables but only first order in the G ones. 

Remark 5.15. Theorem 5.11 is similar in spirit to the reduction theorems for 
classical mechanical systems, like Theorem 2.7. Still there is a noticeable technical 
difference between both types of results. In the continuous case, even in the gen- 
eralized context, the choice of a connection on the principal bundle tt serves the 
dual purpose of constructing a model for the reduced space via the diffeomorphism 
and determines a horizontal / vertical splitting of the variational principle and 
equations of motion. In the discrete context a continuous connection serves the 
same splitting purpose but a discrete connection Ad is used to construct the model 
reduced space. It would be interesting to see if there is any advantage in using two 
different connections in the reduction of continuous mechanical systems. 

Example 5.16. Continuing the analysis of the discrete mechanical system intro- 
duced in Example 5.1, we apply Theorem 5.11 to find the equations of motion of 
the reduced system. The reduced variational constraint is 

A, =d^((rfe,0))(P(,,^o)) = rf'r((rfc,0))((9,|(^^_oj-rfeay|(^^_oj)) 
= (5rL )=T,,,(Q/G), 
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and, since 

{ro,wo,ri) (zVd^ $ ,J ((ro, 0), wq, ''i) £ ^ Wo = (1 - b){rl -rl)/2, 

d 

the reduced kinematic constraints arc 

Vd = {(ro, «;o, ri) e (g/G) x G x (Q/G) : w,, ^ {I ^ h){rl - rl)/2}. 

Next we compute the reduced equations of motion. Notice that iS = {0} for this 
system, so that ~ {0} and condition (18) is trivially satisfied, that is, there are 
no vertical equations. 

In order to find the horizontal equations we compute 

DiLd{{rk,hh),Wk,rk+i) = -m((rfe+i - rfc) + (wfe + &(rfe+i - rl)/2)brk)dxk\f^^^ f^^-^ 
D:>,Ld{{rk-i,hk-i),Wk~i,rk) = m{{rk - ru-i) + (wfe-i + b[rl - rl^^) /2)brk)drk\^^ 
recall from Example 5.6 that 

and, from Example 5.9 the expression of the reduced forces. 

F;[{rk,Wk.,rk+i) = - m(l - b){wk + K^^+i - r^)/2)rfedrfc 

Fj^{rk-i,Wk-i,rk) = m{l - b){wk-i + b{rl - rl_^)/2)rkdrk\^^. 

Then, (p = mUkdrk \^ for 

Uk ■■= - {{rk+i - rk) - {rk - ffc-i)) 

- rk{wk - Wk-i + biirl^, - rl) - {rl - rti))/2), 

and, since V ~ T{Q/G), condition (16) says that Uk — 0, which is the horizon- 
tal equation of motion for the reduced system. Thus, the reduced evolution is 
determined by the system of equations 

(21) C/fc=0 and Wk = {I - b)lrl_,, - rl) /2. 

Using the second equation to eliminate the w dependence in the first expression we 
obtain 

(22) = {{vk+i - rk) - {rk - rk-i)) + rk{{rl+, - rl) - {rl - rl_,))/2. 

From this equation of degree two in rk+i the evolution of this variable is obtained 
and, using the second equation of (21) the dynamics of Wk is determined. 

Remark 5.17. Notice that, depending on the value of the parameter b, the reduced 
system constructed in Example 5.16 is forced or not. Besides making the reduced 
system unforced, the value 6=1 also gives a very simple dynamics to the w. 
variables. This will be a characteristic that we explore in more detail in Section 10 
when we consider systems of Chaplygin type. 

Back in the general setting, conditions (16) and (18) in Theorem 5.11 establish 
the equations of motion of the reduced system. However, they are explicitly written 
in terms of q.. We will see in Section 6 that, in fact, they can be defined in an 
intrinsic manner, in terms of objects defined on the reduced space. 
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6. Intrinsic version of the reduced equations of motion 

In this section wc write the horizontal and vertical equations in terms of the 
reduced system. In fact, the equations of motion will be given as conditions on 

(2) 

morphisms of vector bundles on the second order reduced manifold Qq := G x 
(Q/G) Xq/g G X {Q/G), where the fibered product is taken over the maps p2 : 
G X (Q/G) Q/G and op^-.Gx {Q/G) Q/G. 

It is convenient to consider the space :— QxGx{Q/G)xq/qGx{Q/G) with 

- (2) 

the G-action Ig ° iqo,wo,ri,vi,r2) {lf{qo),lgiwo),ri,vi,r2), so that Qq^ /G = 

(2) 

Qq . Additionally, we define the maps 

Fi:Q)-.'^Q by Fi(go, wq, ri, i^i(go, wq, ri), 



and F2 : Q^^ ^ G by 



^2(go,Wo,ri,wi,r2) /^(^^(^^^^^ ^^^(lii) if vi=p{qi,wi), 

where t : Q Xq/q Q ^ G is defined by T[q, q') = g ii lf{q') = q- It is easy to check 
that Fi and F2 are G-equivariant. 

6.1. Horizontal equations. In order to give an intrinsic meaning to (16), we 
consider the following commutative diagram, 



(23) R^J^p*T{Q/G) T{Q/G) 




where 4> is defined by 

0(go, WQ,ri,vi,r2,6ri) ■.={DiLd{qi,wi,r2) o h"^^ + DzLd{qa,wa,ri) 

+ F-ivi,r2) + F+{piqo,wo),ri)){Sri) 

for qi := Fi{qo,wo,ri) and wi := F2{qo,wa,ri,vi,r2). Also, G acts on plT{Q/G) 
hy lg^^^'^^'^\qo,wo,ri,vi,r2,Sri) := {l'^{qo),l^{wo),ri,vi,r2,dri) and trivially on 
R. 

Lemma 6.1. The map <j} is a G-equivariant morphism of vector bundles. 

Proof. From the explicit definition, (p is a morphism of vector bundles. 

" (2) 

In order to check the G-equi variance, for {qQ,wo,ri,vi,r2) S Qg > we let qi := 
Fiiqo,WQ,ri) and wi := F2iqQ,wo,ri,vi,r2). Then, writing 

(j)i{qo,Wo,ri,vi,r2,Sri) := {DiLdiqi,Wi,r2) o h'^^ + D3Ld{qQ,Wo,ri)){Sri) 

(t)2{qo,Wo,ri,vi,r2,5ri) {F~{vi,r2) + F^{p{qo,wo),ri)){dri), 

we have (j) — <j}i+ (j)2. From the definition, it is immediate that 02 is G-equivariant, 
so we concentrate on (j)i . 

'^i('s'^^'^^'^^(9o,wo,ri,wi,r2,(5ri)) {qo),!^ {wQ),ri,vi,r2,5ri) 

(24) 

= i^i id (/,^ (<?i ) , ( ) , r2 ) o /i',' (^ r 1 ) + i?3^d (/,^ (<?o ) , (^0 ) , ri ) (<5r 1 ) , 
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where we used the G-equivariance of Fi and F2 in the first term of the last identity. 
By Lemma 5.2, we have DiLd{lf {qi) , if {wi) , r2)l'^'^ = DiLd{qi,wi,r2), so that 

D,UlQiq,),lfiw,),r2)oh'?^''^^ = D,UlQ{q,)jf{w,),r2)ol'^Qoh'^^ 

= DiLd{qi,Wi,r2) o h'^K 

Also, from Ld{lf {qo),lf {wo),ri) = Ld{qo,wo,ri), we obtain 

D3Ld{l'^{qo),lg{wQ),ri) = D3Ld{qo,Wo,ri). 

Replacing the last two identities back in (24), the G-equivariance of (pi follows. □ 

By Lemma 6.1, (j> defines a morphism of vector bundles cj) : {p%T{Q/G))/G R. 
Since _(p5r(Q/G))/G = plT{QIG), where p2 : Q^?^ ^ Q/G is the projection, we 
have (j) ■ P*2T{QIG) ^ M. Concretely, if {qQ^wa) G p~^(wo), 

cj){vo,ri,vi,r2,5ri) = 0(go, Wq, ri, ui, r2, fri). 

We extend diagram (23) to the following commutative diagram, where T) C 
T{Q/G) is the subbundle introduced in Theorem 5.11. 



T{Q/G) 




Q/G 



Proposition 6.2. Condition (16) in Theorem 5.11 is equivalent to 
(25) 



lp*(X)) 



= 0. 



Proof. For Vk and defined as in Theorem 5.11, the explicit form of coincides 
with that of (j) defined in (15) and the vanishing condition (16) coincides with (25). 



□ 



6.2. Vertical equations. Consider the commutative diagram 



(26) 



Q X 




where 



■4){qo,wo,ri,vi,r2,^i) := {D2Ld{qo,wo,ri)wQ ^ - D2Ld{qi,wi,r2)wi "^)(Ci) 
for qi ;= Fi{qo,wo,ri) and wi := i^2('?o, ^"0, ''i, ^i, ?'2)- 
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" (2) 

Notice that F^Q x q) = Qq x g. In addition, G acts on M trivially and on 
Qg) X by ;3^^""^<Zo,^^o,ri,z;i,r2,a) ^^(go), /g^K), «i, r2, 
Lemma 6.3. The map ip is a G-equivariant morphism of vector bundles. 
Proof. From the explicit definition, -0 is a morphism of vector bundles. 

" (2) 

In order to check the G-equivariance, for {qo,wo,ri,vi,r2) £ Qg , we let qi :— 
Fi{qo,WQ,ri) and wi -.^ F2{qQ,wo,ri,vi,r2). Since 

ip{qo,wo,ri,vi,r2,^i) = D2Ld{qo,Wo,ri){^iWo) - D2Ld{qi,Wi,r2)i^iWi), 

using the G-equivariance of Fi and F2, we have 

ipi^g^ ''°iqo,wo,ri,vi,r2,^i)) =-;/'(/^((jo), ^^(^o), ri, wi, rj, 

(27) =D2Ul^ (go ) , If (wo ) , ri ) (/™ i^i)lf (wo)) 

-D2Ld{lf{qi),lf{w,),r2){l^''mf{w,)). 
On one hand, recalling the notation of Remark 5.10, we have 

(28) lg'^iii)lf{wo) = g^ig-^gwog-^ = 56^^0.9"' = dlf{wo){^iWo). 

On the other, the G-equivariance of D2Ld proved in Lemma 5.2 together with (28), 
applied to (27) lead to the G-equivariance of ip. □ 

— ^ /2) 

By Lemma 6.3, '0 defines a morphism of vector bundles ip : {Qq x g)/G K 
over Q^Q . Exphcitly, for any [q^^wo) € p~^{vo) and ^1 S g, 

(29) _ 

■ipi.[{qo,Wo,ri,Vi,r2,^i)]) ='(/'((7o, Wq, ri, wi, ra, fi) 

^{F)2Ld{qa, Wq, ri)wQ^ - D2Ld{qi,Wi,r2)w^^) (^1), 

where qi := Fi{qo,wo,ri) and wi := F2{qo,wo,ri,vi,r2). 

A simple computation shows that, because T> C TQ is G-invariant, F*g^ C 
Qq^ X g is G-invariant. Therefore {F^q'^)/G C {Qq'' x g)/G is a vector subbundle. 

We collect the different objects in the following commutative diagram, which is 
the quotient of (part of) (26) 



-,(2) 




g)/G- ^(^^rg^)/G 



Proposition 6.4. Condition (18) in Theorem 5.11 is equivalent to 
(30) |(j..g.)/G = 0. 

Proof. For Vk and rk defined as in Theorem 5.11, the explicit form of "0 coincides 
with that of ij} defined in (17) and the vanishing condition (18) coincides with (30). 

□ 

Corollary 6.5. Any one of the four equivalent conditions of Theorem 5.11 is equiv- 
alent to (ufe,rfe-|_i) G Vd and conditions (25) and (30) are met for all k. 
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7. Reconstruction 

Given a discrete curve q. in Q and its image (u.,r.) in G x [Q/G), Theorem 5.11 
establishes an equivalence between q. being a trajectory of the original system (i.e., 
a solution of the original dynamics) and {v.,r.) being a solution of the reduced 
dynamics. In this section we study the reconstruction problem, that is, given 
a curve {v.^r.) in G x (Q/G) that satisfies adequate conditions (constraints and 
equations of motion), is it possible to find a trajectory q. of the original system 
that projects to {v.,r.)l 

Consider the following construction. Given a discrete curve {v.,r.) £ G x (Q/G) 
and qk G Q (one value of k) such that TT{qk) = Vk ~ p'^^'^{vk), if Wfe = p{qk,'Wk) we 
define 

(31) Uk ^?{qf,.qk)i'^k) e G and qk+i Fi{qk,Uk,rk+i). 

^^""'^ ^riqk.ifiq,)^^f(^~k)) ^?iq.,qk),-^ilfi^'')) = the previous 

construction is independent of the chosen representatives of Vk- Therefore, given 
qo & Q with 7r((7o) = t'oi applying the construction iteratively defines a unique 
discrete curve q. in Q. The following result establishes the properties of q.. 

Theorem 7.1. Consider {qo,qi) G T^d o,nd {v.,r.) e G x {Q/G) such that 7r(go) = 
To, 7r(gi) = ri, vq = p{qQ,Ad{qa,qi))- If {v.,r.) satisfies {vk,rk+i) € Vd and 
conditions (25) and (30) for all k, then the discrete curve q. constructed by (31) 
from qo is a trajectory of the system on Q whose image by T is the curve {v.,r.) 
and satisfies qo = go and qi ^ qi. 

Proof. The curve q. is a lifting of the curve (w.,r.) to Q; that is, = T^iqk) 
and Vk ^ p{qk,Ad{qk,qk+i))- Indeed, from (31), 7r(gfc) = 7r(/g^_^ (/i*"' (rfc))) = 
7r(/i^'°"^ (rfc)) = Tfe. On the other hand, 

Ad{qk,qk+i) ^ Ad{qk,l2^{hf {rk+i))) = UkAd [qk, hf (rk+i)) =Uk, 

so that 

p{qk,Ad{qk-,qk+l)) = p{qk,Uk) = PiQk,lfl^q^^q^-){Wk)) 

= Pi''?(q,,q,)ilk)^^k) = p{qk,Wk) = Vk, 

which completes the argument, that is T{qk,qk+i) = {vk,rk+i) for all k. 

The lifted curve satisfies the first initial condition because, by construction, rg 
is lifted to qo- Since vq = p{qo,Ad{qo,Qi)), according to (31), uq = Adiqa^Qi) so 
that ri is lifted to qi = l^^^g^ qi)i^'d' (''"'i-))- Using (8), we conclude that qi = qi. 
Next, we see that {qk, qk+i) & T^d for all k. By the G-invariance of Vd, 
{qk,qk+l) e'Dd <^n{qk,qk+i) eVd/G ^ (^(gfc, Qfe+i)) G (2?d/G) 
<^ {vk,rk+i) G Vd, 

which holds by hypothesis for all k. 

The only thing left to do is to check that q. is a trajectory of the discrete me- 
chanical system. By hypothesis, {v., r.) satisfies conditions (25) and (30). Then, by 
Corollary 6.5 conditions (16) and (18) hold. But, since the relationship among q., 
v., r. and w. is precisely that of the statement of Theorem 5.11 and we proved that 
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condition 4 holds, we conclude that q. satisfies condition 1 in Theorem 5.11, hence 
it is a trajectory of the system on Q. □ 

Example 7.2. The last step to complete our analysis of the system {Q, Ld, T), T>d) 
introduced in Example 5.1 is to consider the reconstruction of the evolution of the 
original system given a trajectory (r.,w.) of the reduced system compatible with 
some initial data {qo, qi) G Vd- 

According to Theorem 7.1 we use (31) to construct the trajectory q.. Since 
Vk = (rfe,Wfc) = p((rfc,0),w;fc) and r((a;fc, i/fe), (r^, 0)) = yu (remember that Xk = 
TT{xk,yk) = rk), 

'^k = ^Vki'^k) = Wk, 

qk+i = {xk+uyk+i)^l2,{h''d'"''"\-rk+i)) = (rk+i,yk + b{rl^, - xl)/2 + Wk) 

that expresses q. in terms of the known data (r. , w.). Simplifying we obtain 

ixk,yk) = irk,yo + (rl- rl)/2) for all k, 

with rQ = xq, ri = xi. Notice that the resulting trajectory is independent of the 
parameter b chosen to do the reduction, as it should be. 

Remark 7.3. The reduction and reconstruction techniques developed so far can 
be applied in the case where the configuration space is the symmetry group (acting 
by left multiplication), a problem that has already been studied in [12, 1-3, 2'd]. In 
this context, Theorems 5.11 and 7.1 allow us to re derive Theorem 3 and Corollary 
4 of [23]. 

8. NONHOLONOMIC DISCRETE MOMENTUM 

For discrete or continuous holonomic mechanical systems the presence of con- 
tinuous symmetries automatically leads to the existence of conserved quantities 
(momenta), due to Noether's Theorem. In the nonholonomic case, this is no longer 
true, essentially due to the behavior of constraint forces. Instead, in this case, one 
obtains an equation that describes the evolution of momenta over the trajectory of 
the system. Below we discuss the relationship of the discrete momentum evolution 
equation and the dynamics of a discrete mechanical system with symmetries. 

Definition 8.1. Given {Q,Ld^'D,'Dd) with symmetry group G, the nonholonomic 
discrete momentum map is the application Jd Q y- Q ^ (fl^)* defined by 

(32) Jd{qo,qi){qo,0 ■= -DiLd{qo,qi)£,Q{qo). 
for aU ((70,0 e . 

Given any section ^ G T{q'^), define the map [Jd)^ : Q ><■ Q ^M.hy 

{Jd)^{qo,qi) ■■= Jd{qo,qi){qo,£.{qo))- 

Straightforward computations give the following result. 
Lemma 8.2. For ^ G g and {qo, qi) ^ Q x Q, 

(33) DiLd{qo,qi){£,Q{qo)) = -D2Ld{qo, qi){^Q{qi))- 
Also, if{qo,Wo,ri) = $^^(90,91) = {qo,Ad{qo,qi),TT{qi)), then 

(34) - DiLd{qo,qi){£,Q{qo)) = {D2Ld{qo,Wo,ri)wQ^){£,). 
In particular, Jdiqo,qi) = D2Ldiqo,Wo,ri)wQ^ . 
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Remark 8.3. By definition, when Q x Q and q are seen as bundles over Q with 
respect to the projection on the first variable, Jd is a bundle mapping (that is, it 
maps fibers to fibers). J. Cortes defines J^^ in [11, pp. 154-5] in a slightly difi'erent 
way: J2^iqQ,qi){qi,0 -D2-^d('7o, '7i)'?q(9i), which is a bundle map when the 
projection on the second variable is considered in Q x Q. In any case, due to (33), 
both maps are, essentially, the same. Also (Jd)^ = i^d^)^' f^^' ? ^ r(0'^). 

A trajectory of a discrete mechanical system is determined by the Discrete 
Lagrange-D'Alembert Principle (Definition 3.3). When the variational constraints 
are decomposed by 2? = 5©H, it is possible to decompose all admissible variations 
into horizontal and vertical variations, in the sense that they belong to H or S. It 
is also possible to decompose the variational principle accordingly. The following 
result, whose proof is obvious, makes a precise statement. 

Proposition 8.4. Let q. be a discrete curve in Q. Then, the following conditions 
are equivalent. 

(1) q. satisfies the variational part of the discrete Lagrange-D'Alembert princi- 
ple (Definition 3.3). 

(2) q. satisfies 

dSd{q.)i6q^) = and dSdiq.){Sq'^) ^ 

for all pairs of variations vanishing at the end points dq'^ and 5q^ in S and 
% respectively. 

The next result relates the vertical variational principle, the evolution of the 
nonholonomic discrete momentum and condition (30). 

Theorem 8.5. In the context of the theorem 5.11, if q. is a discrete curve in Q, let 
rk = '^{qk),Wk = Ad{qk,qk+i) and Vk = p{qk,Wk)- Then, the following conditions 
are equivalent. 

(1) q. satisfies the vertical variational principle. That is, dSd{q ){5q^) = for 
all vanishing end point variation 5q^ S Sq^ for all k. 

(2) Condition (30) is satisfied. 

(3) For all sections ^ € T{g^), 

(35) {Jd)^{qk,qk+i)-{Jd)^{qk-i,qk) = -L'i£d('7fe-i, <?/=)(!(%) -|('7fc-i))Q('7fc-i)- 

Proof. • 1 <^ 2. This is Proposition 6.4 and Lemma 5.12. 

• 2^3. For ^ e r{Q^), using the definition of (Jd)|, equation (35) is 
equivalent to 

(36) - DiLd{qk,qk+i){i{qk))Q{qk) = -~DiLd{qk-i,qk){i{qk))Q{qk-i)- 

Using now (34) with ^ = C('Zfe) a-nd recalling the explicit formula (29) for 
we obtain (30). 

□ 

Remark 8.6. During the proof of Theorem 8.5, we saw that the discrete nonholo- 
nomic momentum evolution equation was equivalent to (36), that is simpler to use 
in practice and, by the same Theorem is equivalent to any of the other conditions 
in the statement. 
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Remark 8.7. The discrete nonholoiiomic evolution equation (35) was first obtained 
by Cortes and Martinez in [10, Thm. 5.3], where they prove that any solution 
of the discrete Lagrange-D'Alembert equation satisfies (35). In the context of 
groupoids, Iglesias et al. obtain the evolution equation of the discrete nonholonomic 
momentum map in [ l o, Thm. 3.20] 

9. Reduced equations of motion: trivial bundle case 

In this section we consider the case where Q := i? x G, for a manifold R and a 
Lie group G that acts on Q by lf{r, h) := (r, gh). In this case, it : Q ^ Q/G is the 
trivial principal bundle pi : R x G R. The goal is to give an explicit description 
of the reduced system as well as the corresponding equations of motion. 

Let a : G ^ R X G he given by a(/9((ro, /iq), wo)) •= {'''Ot ^^o^WotiQ), then a is 
a diffeomorphism with inverse /3(ro,i9o) := p((ro, e), i^o)- We have the following 
commutative diagram 

QxGx (Q/G) RxGxGxR 

jixid 

pxid xid 

G X {Q/G) RxGx R 

f)xid 

where p*(ro,/io,'?o) (?'o, V^'^o^o), a((ro, ft-o), ^^o) := (j^c^c^^o) and /3 a'^. In 
addition, a useful ingredient is the section s : G QxG given by s(p((ro, ft-o); i^o)) ■= 
((ro,e),/io ^7?o^o)- 

In the current context, by (7), Ad can be written as Ad{{rQ,ho),{ri,hi)) = 
/ii.4^(ro,ri)/io \ where y4.^(ro,ri) := ^^((ro, e), (ri, e)). 

Using a the reduced constraint manifold is 2?^ := (a x id){'Dd), so that the 
reduced kinematic constraint condition becomes {ro,'do,ri) G 2?^ if and only if 
{{ro,e),{ri,doA';j{ro,ri)-^)) G Vd- If Vd is described by equations (j>b{qQ,qi) = 
for all b, the reduced constraint condition becomes 0h((ro, e), (ri, i9o^d(''Oj ''i)^^)) = 
for all b. 

We also have the induced lagrangians := Lrfo (/3 x id) and := Ldo{(ix id), 
so that Ll{ro,iSo,ri) = Ld{{ro,e),'&o,ri). 

" (2) (2) 

Next we characterize the second order reduced manifolds Qq and Qq , that are 
needed to set the equations of motion, according to Corollary 6.5. It is clear that 

'^RxGxGxRxGxR, 

with ((ro,/io), wo,?'i,p((fi,'ii), wi),r2) >-> {rn,h(3,wn,ri,h^'^wihi,r2), and 

Q^i^^ ^RxGxRxGxR, 

with{p{{ro,ho),wo),ri,p({ri,hi),wi),r2) i-^ (ro, ,hQ^Woho,ri,h^^wihi,r2). Next, 

(2) . ~ ~ 

we have to characterize the bundles over Qq where the morphisms </> and ip are 
defined. Instead, it is easier to notice that the section s provides a diffeomorphism 

f21 f2l 

of Qq with its image, so that we can view Qq' as Rx {e} xGxRxGxR. The 
advantage of this approach is that instead of having to work with quotient bundles 

" (2) 

we have to work with the restriction of bundles on Qq to the image of s. 
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We start by obtaining trivialized versions of the maps Fi, Fi and F2. They are 



Remark 9.1. The case when the symmetry group is abelian has some speciahy 
nice features. For one thing, even when tt : Q — Q/G is not trivial, there is a 
diffeomorphism a : G — > (Q/G) x G that is given by a{p{qo,wo)) ■= {TT{qo),wo). 
When, in addition, tt : Q — > Q/G is trivial, this diffeomorphism coincides with the 
one introduced at the beginning of this section. 

— (2) 

9.1. Horizontal equations. Here we write down the mor phism (f> : (g^^ xp^TR) -> 

(2) 

M. or, rather, its realization under the identification of Qq with the image of the 
section s in Qq^- Let 0* be the pullback of (f) to RxGxGxRxGxR. Explicitly, 

0*(ro, /lo, -ffo, ri,di,r2, Sri) := (^((rg, /iq), i9o, ri, p{{ri, e),'di), r2,Sri) 

= (DiLdiiri,hi),wi,r2) o h'^^ + D3Ldiiro,ho),-do,ri) 

+ (F,-)*(ri,^i,r2) + {F+nro,i}o,n)){Sn), 

where hi := ?9o/io(^d(»'o, qi ■= {ri.hi), wi := i^|(ro, /iq, ^^o, J'l, ^^i, ^'2), and 

{F^Y{ri,di,r2) := F^{p{{ri,e), i^i), ra). 

From Definition 5.8, we see that F'^ : p\TR®p'^TR — > R is a morphism of bundles. 
Explicitly, 

Fd{ro,-dQ,ri){5rQ,5ri) -.^ 

D2Ldi{ro, e), 1^0, ri)dAd{{ro, e), {n, hi)){h^-^>''\6ro), h^'-^'^'HSn)), 

with hi as above. Even more, we can write 

dAd{{ro,ho),{ri,hi)){6ro,6ho,Sri,Shi) := hiDiAd{ro,ri){6r„)h^^ 

- hiA'diro, ri)h^^{5ho)h-^ + hiD2A'd{ro, ri){6ri)h-^ + ShiA'diro,ri)h^\ 

and we notice that this expression is written only in terms of 

We can also write the horizontal lift explicitly. In order to do so, we view the 
vector bundle W© H C TQ as the graph of a bundle map M : plTR -5- p^TG over 
Q = RxG. Therefore, 



Fi*(ro,/io,i?o,r-i) := Fi((ro, /iq), i9o, = {ri,i}ohoMro,ri)y^), 
Fl{ro,ho,'0o,ri,di,r2) := Fi((ro, /iq), ^?o, ^i, p((ri, e), 7-2) 
= {ri,i)oh,{AUro:ri))-'), 



and 



F2{ro,ho,'dQ,ri,di,r2) := -F2((''o, /io), 1^0, '"i, e), ^2) 

= ^ohoiA'diro, ri)y')^iA'd{ro, ri) V'l^o 



(<5r) := {Sr,M{r,h){5r)) e T^r.h)Q- 



Finally, we have 

4)*{ro,-do,ri,di,r2,Sri) := (j)\rQ,e,'do,ri,-di,r2,5ri) 



(37) 




+ D2L'^{ri,§i,r2) o {diM{ri,e) - M{ri,e)di) 
+ (/'-)*(ri,i?i,r2) + (F+)*(ro,i^o,ri))(<5ri). 
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On the other hand, we have that the representation of P2'D over RxGx RxGx R 
is the vector bundle P* := p^T^- 
Finally, condition (25) becomes 

(38) 0*(ro,7?o,ri,i?i,r2,^ri) =0 

for aU 6ri such that (fri) g P(r,.e). 

Alternatively, if (wi, . . . ,u!k) = 1^°, condition (38) becomes 

0*(ro,i9o,?-i,i?i,r-2,5ri) = Aatj''(ri, e)((5ri, M(ri, e)((5ri)) 

^ — ^a — 1 

for all dri £ Tr-^ R and where Aa £ K are unknown. 

9.2. Vertical equations. Here we write down the morphism ■ {Qg ^Q)/^ ~^ ^ 

(2) 

or, rather, its realization under the identification of Qq with the image of the 
section s in Let i/j* be the puUback oi ip to RxGxGxRxGxR. Explicitly, 

ip*{ro,ho,{)o,ri,di,r2,^i) := ■(/'((ro, /iq), i^o, ^i, p((ri, e), i?i), r2, ^i) 

= {D2Ld{iro,ho),'do,ri)wQ^ - D2Ld{{ri,hi),wi,)wY'^){(,i), 
where hi := ?9o^o(-4^(?'o, ?'i))"^ and wi := F^{ro, ho,do,ri,'di,r2). 

Since, V* := V'*|i?x{e}xGx_RxGx_R, 

■ip*{ro,-do,ri,^i,r2,^i) = ?^*(ro, e, ??o, ''i, i?i, ?'2, 6) 

(39) = (i?2id((ro,e),i?o,r-i)i9^' -7^2id((ri,/ii),wi,r2)wr')(6) 
= {D2L*dirn,^n,n)^^^ ~ hiD2L*diri,^i,r2)^^^h^^)i^i), 

where 

(40) /ii := i?o(^d(?'o,ri))^i and wi := F^{ra,e,i}o,ri,i}i,r2). 

The other required ingredient is the pullback of F^g^ , that is, {FI)*q^ . Then, 

(Ff)*fl^ = {{ro,ho,^o,n,^i,r2,Ci) e Q^'' >< ■ ((^i, /ii), (0, 6'^)) e T^ir.M)} 

= {(ro,/io,i?o,?'i,i?i,r2,6) G <9g^ x g : ((ri, e), (0, /ij^^^i/ii)) g 2?(ri,e)}, 
where /ii is as in (40). Therefore, condition (30), becomes 

(41) V^*(ro,i?o,n,i9i,r2,a) -0 

for aU ^1 g such that (0, h^^S,ihi) g where /ii is defined by (40). 

Alternatively, if (wi, . . . ,ujk) — T^°, condition (41) becomes 

(42) 7/;*(ro,z?o,ri,7?i,r2,Ci) - V"" A,L.^(ri,e)(0,/ir'a/ii) for all S fl, 

^ — ^a—l 

where hi is defined by (40) and A^ g K are unknown. 

9.3. Reconstruction. Given initial conditions on Q x Q, the reconstruction pro- 
cess in the current case goes as follows. Let go = i^o-,ho) and 51 ~ (fi,/ii) with 
{qo,Qi) g I'd, and a discrete curve (r., I?.) g i? x G. If (r. , i?.) satisfies the constraint 
P^, the equations (38) and (41), and ro = ro, ri ~ fi and i9o = ^o^^'ii-^dl^o, ^1), 
by Theorem 7.1, the lifted curve q. constructed from go using (31) is a trajectory 
of the original mechanical system on Q x Q with qo = qo and qi = qi- In the trivial 
bundle case, if qk = {rk,hk), 

(43) qk+i = {rk+i,hk'dkAd{rk,rk+i)~'^) for aU /c. 
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9.4. An example of reduction in the trivial bundle case. In this section we 
specialize the previous discussion to the discrete mechanical system {Q, Ld,'D,'Dd) 
where Q := x x with coordinates (x, y, 9, 0) — we consider := ]R/27rZ, 
and identify operations in with regular addition in M} — . 

V, -.^{d^^de + Acosiq'i')d^, + Asin{q^)dy) = {uj,{q),i02{q))° C T^Q, 

where uji{q) := dx — r cos{q'^)d9 and uj2{q) := dy — r sm{q'^)d9, and 

Vd {(go, gi) eQxQ-.q^-q^,^ A{ql - q^oKcos^qt + qf )/2) and 

qf - ql = A{ql - ql){sH{qt + gf )/2)}, 

for TO, A, I and J positive constants. This system can be obtained as a discretization 
of the classical mechanical system formed by a vertical disk of radius A, mass to 
with inertia momenta / and J, rolling without slipping on a horizontal plane. 
We consider the Lie group G := x acting on Q by 

The corresponding lifted action is lj'^{q,v) = {l'^{q),v). We are in the trivial 
bundle case analyzed above because tt : Q — > Q/G is the trivial principal bundle 
P2 '■ G X ^ with structure group G, which is a symmetry group of the system. 

We consider various subbundles of TQ: the vertical bundle = {dx,dy,de), 
the intersection bundle iSg = {de+Acos{q'^)dx + Asm{q'^)dy) and the corresponding 
complements Hq := {d^), Uq := {dx,dy) and W := {0}. Following Definition 2.5, 
the previous bundles induce a nonholonomic connection A whose horizontal space 
is HorA ■= H- The corresponding horizontal lift h"^ : TrS^ — Hq C TqQ for 
q G 7r~^(r) is h'^{rdr) = "rS^; notice that, in this case, the bundle map M that 
describes the horizontal bundle as a graph is identically 0. 

By Proposition 4.12, since 

Ad{q^,q,) = {ql,qlql){q^,qlql)-' = [q^ - q^^ql - qlqf - q^) £ G 

satisfies (7), there is a discrete connection with level 7(5) = e whose discrete connec- 
tion 1-form is Ad- Notice that our choice of Ad has ^^(rg, ri) = e for all ro, ri G . 
The corresponding discrete horizontal hft is h^^'iri) = {qo, (qg , , qg, ^i))- 

The reduced space in this case is G x [Q/G) ~ S'^ x G x S*^ and the reduced 
second order manifold is Q^q ^ x G x x G x . The trivialized reduced 
lagrangian is 

L'dMo,n) = + i^yf) + ^-Kf + rof. 

Next we describe the horizontal and vertical equations of the reduced system, 
but first we notice that since Ad only depends on the group variables, it turns 
out that dAd annihilates horizontal vectors and the mixed curvature Bm vanishes. 
Consequently, there are no reduced discrete forces on the system. 

From (37), we derive 

0*(rfc_i,i?fc_i,rfc,79fc,rfc+i,5rfc) = (J((rfe - rk-i) - {rk+i - rk))dr\r^){5rk), 
and, noticing that W — {0} implies that & = TS^, the horizontal equation is 
(44) (rfc - Tk-i) - (rfc+i - rfc) = 0. 
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From (39), we obtain 

V;*(rfc_i,i9fc_i,rfc,z?fc,rfc+i,Cfc) = {m{dl_^ ~ dl)dx\,. + m{-&l_^ ~ ^l)dy\ey 

+ I{K-i~K)deio){ik), 

where e = (e^, e**, e^) is the identity element in G. Since the right hand side of (42) 
is Ai(da;|e^ — Acos{riS)d6\^e) + \2{dy\ey — As\n(ri;)dO\^e), we obtain the following 
vertical equations 

m(^Li-^|) = A2, 

liK^i - O - - Ai^cos(rfe) - A2Asin(rfe). 
On the other hand, the reduced kinematic constraint equations are 

^fc=^^feC0s((r,+r,+i)/2), 
i?I=A<.sin((r,+rfe+i)/2). 
Next, we find the reduced dynamics by solving the reduced equations. From (44) 

(46) Tk = (ri - ro)fc + tq for all fc £ N U {0}. 
From the vertical equations we obtain 

Plugging (45) into this last equation (for fc — 1 and k) and simplifying we get 

^ I + mA^ cos{{rk-rk-i)/2) ^g 
^ / + 771^2 cos((rfc+i -rfc)/2) 

Using the horizontal equation (44) we conclude that = t^q for all fc. Thus, 
7?^=A<cos(^(ri~ro)fc+^^ii^^ and = A< sin (^(n -ro)fc+ ^^^^^^ 

This last expression completes the description of the reduced system's dynamics. 

The reconstruction of the trajectories of the original system in Q x Q is done as 
described in Section 9.3. By (43), the reconstructed trajectory q. is 

(47) qu+, = {{qt+,.ql^„qU,),qt+,) = {{ql , q^ ql)d u , r , 
SO that {ql+i,ql+i,qi+i) = (^o , 9o) 11^=0 '^i or, 

,sin((gf - q^)k + qf,) - s\n{qf,) 



ql^, = gj? + A{ql ql) cos((gf - ql)/2)- 



sin(gf - q^) 



- 4 + Aiq^, - cos((,f - ^^^y, f-<lt)~cosiiqt-qt)k + r.) 



sm{qf - ql) 



ql+i = iql-ql)k + ql. 

Also, from (47) and (46), 

qU, = {qt-qt){k + l) + qt- 
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10. Reduced equations of motion: Chaplygin case 

In this section we specialize Theorems 5.11 and 7.1 to the case where the original 
system is a discrete mechanical system with Chaplygin type symmetries. In this 
case we go beyond the result of equivalence between the discrete mechanical system 
in Q and a dynamical system in Gx{Q/G) to obtain an equivalence between discrete 
mechanical systems on Q and Q/G. 

Definition 10.1. A symmetry group G of [Q, Ld,'D,'Dci) is a Chaplygin type sym- 
metry group if it satisfies the following conditions. 

(1) TQ = V^®V and 

(2) Vd defines an affine discrete connection Ad on the principal bundle tt : Q ^ 
Q/G. 

Notice that, by definition, the condition TQ = V'~^ © I? is equivalent to the fact 
that the decomposition (2) has the form 

(48) TQ = {0} ©J^© {0} 

As in the general case, this decomposition defines a connection ^ on Q Q/G, 
whose horizontal space is "H = P. Condition 2 in Definition 10.1 requires that, for 
every {qo,qi) there exists a unique g €G such that (90,^^-1(91)) G 'Dd- 

Example 10.2. G := K is a Chaplygin type symmetry group of the discrete 
mechanical system {Q, Ldj'Dj'Dd) defined in Example 5.1. The first condition is 
clear from the decomposition of TQ that appears in Example 5.6; the second point 
corresponds to the discrete connection Ad considered in Example 4.17 when 6=1. 

10.1. An inclusion. Let y : {Q/G) x (Q/G) ^ G x (Q/G) be defined by 

y{r(i,ri) := {p{qQ, e),ri), 
where qq S 7r~^(ro) and e is the identity of G. 

Lemma 10.3. The application y is well defined. Even more, if s is a local section 
ofir-.Q^ Q/G, then y{ro,ri) = (p(s(ro), e), ri). 

Proof, y does not depend on the choice of go G '"'"^(^'o) by the G-invariance of p. 
The expression of y in terms of s holds because s(ro) € 7r~^(ro). This formula also 
shows that 3^ is a smooth map. □ 

Using y we transport the existing structure on G x (Q/G) to {Q/G) x (Q/G). 
More precisely, we define the (forced, unconstrained) discrete mechanical system 
(Q/G, Ld, Fd) where Ld := y*{Ld) = Ld o y and the discrete force Fd := y*{Fd). 
We also define 2? ~ P as a subbundle of T{Q/G), but notice that, V = T{Q/G) 
in the Chaplygin case, due to (48). 

Lemma 10.4. If A is a connection on the principal bundle tt : Q — ?> Q/G and 

go e 7r"i(ro), then 

(49) dLd{ro, ri){Sro,Sri) = D^qo, e, ri){h'i" {Sro)) + DMqo,^, ri)(5ri). 
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Proof. If s is a local section of Q — s> Q/G with go = s(?'o) then Ld{ro,ri) = 
Ld(p(s(ro),e),ri) = Ld{s{ro),e,ri), so that 

dLdiro,ri){Sro,Sri) ^ {dLd){s{ro),e,ri){ds{ra){5ro),0,Sri) 

(50) 

^DiLd{s{ro),e, ri){ds{ro){Sro)) + D3Ld{s{ro),e, ri){dri). 
Since tt o s = idq/Q, we have 

5ra = (i7r(s(ro))ds(ro)((5ro) = d7r(s(ro))i?or^(ds(ro)(5ro)), 

thus iJor_4((is(ro)(5ro)) = (Jro). Also, as DiLd{s{ro),e,ri) vanishes on ver- 
tical vectors by the G-invariance of Ld, 

DiLdis(ro),e, ri){ds{ro){Sro)) =L)iid(s(ro), e, ri)(iJor^((is(ro)((5ro))) 

=D,Ld{s{ro),e,n){h'^-'>)^Sro)) 

which, replaced in (50) leads to (49). □ 

Lemma 10.5. Let (v.,r.) be a discrete curve in G x (Q/G) such that = p{qk,e) 
for some qk £ Q and all k. Then 

dSd{v.,r.){Sv.,Sr.) ^ dSd{r.){Sr.) + (^d"(^fc,rfc+i) + F+(rfc_i,rfc))(frfc) 

fc=i 

for all vanishing end points variations Sr. and 

5vk = dp{qk,e){h''^ {Sr,,),dAd{qk,qk+lW {Srk),h'^''^' (Srk+i))) 

for all k. 

Proof. By definition, using the explicit form of dv, 

dSd{v.,r.){bv.,Sr) =Y.k^o {DiLd{qk,e,rk+iW'' (Sr^)) 

+ D2Ld{qk,e,rk+i){dAd{qk,qk+i){h'"'[5rk),h'"'+^[5rk+i))) 
+ DsLdiqk, e,rk+i){Srk+i)) . 
Using Lemma 10.4 and the decomposition dAd ~ DiAd + D2Ad we obtain 
dSdiv.,r.){6v.,Sr.) = dS{r.){6r.) 

+ Z^fe=o i^d {qk,e,rk+i){Srk) + FJ{qk,e,rk+i){6rk+i)), 
and the result follows. □ 

Lemma 10.6. Let {v.,r.) be a discrete curve in G x (Q/G) such that = p{qk,e) 
for some qk G Q and all k. Then (w.,r.) satisfies condition (25) if and only if 

(51) DiLd{rk,rk+i) + D2Ld{rk-i,rk) + F+{rk-i,rk) + Fj;{rk,rk+i) G (PrJ° 
for all k. 

Proof. Considering the form of Vk, the vanishing condition (25) immediately trans- 
lates, via y to the vanishing condition (51). □ 
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10.2. Reduced dynamics. 

Theorem 10.7. Let G be a Chaplygin symmetry group of (Q, Ld,'D,'Dd), q- be 
a discrete curve in Q, and rk := T^{qk) be the corresponding curve in Q/G. If 
{qk,qk+i) G T^d for all k, the following statements are equivalent. 

(1) q. satisfies the variational principle dSd{q){Sq.) = for all vanishing end 
points variations 5q. such that Sqk G for all k. 

(2) q. satisfies the Lagrange-D'Alembert equations (3). 

(3) r. satisfies the variational principle 

dSd{r-){Sr.) = - ^(F,-(rfc,rfc+i) + F+(rfc_i,r,))(<5rfc) 
fe=i 

for all vanishing end points variations Sr. with Srk G Tr^{Q/G) and where 

N-l 

Sdir.) := ^ Ld{rk,rk+i). 

k=0 

(4) r. satisfies the following equation, for all k. 

DiLd{rk,rk+i) + D2Ld{rk-i,rk) = -(i^^+(rfc_i, r^) + F^^ {rk,rk+i)). 

Proof. Since S = {0}, by Theorem 5.11, we see that each one of the pomts 1 and 2 
of the present result is equivalent to any one of 

(i) Item 3 in Theorem 5.11 holds for all horizontal variations, {£,k = for all 
k). 

(ii) For all k {vk,rk+i) S Vd and condition (16) holds. 

As HorAa = 'Dd, {qk,qk+i) e Vd if and only if Wk = Ad{qk,qk+i) = e. Then, 

{vk,rk+i) eVd^ {p{qk,e),rk+i) £ $^,(23d/G) ^ {qk,h^ {rk+i)) G Pd, 

which always holds by definition of Ad- Therefore 2?^ = G x (Q/G), so that we can 
drop the reduced kinematic constraint condition from i and ii. 

By Lemma 10.5 the variational principle that appears in i is equivalent to the 
one in point 3 of the statement. Similarly, as V — T{Q/G), Lemma 10.6 and 
Proposition 6.2 show that ii is equivalent to the one in point 4 of the statement. □ 

Remark 10.8. Notice that the reduction of a discrete mechanical system with 
Chaplygin type symmetry results in an unconstrained discrete mechanical system 
but with external forces given by in the previous theorem. A similar analysis 
has been done in the more general groupoid setting in [15]. 

Example 10.9. As we noted above, G = M is a Chaplygin type symmetry group 
of the discrete mechanical system introduced in Example 5.1. We notice that the 
construction described in this section corresponds to performing the reduction us- 
ing the discrete connection from Example 4.17. Furthermore, the inclusion 
y from section 10.1 was essentially already present in the analysis of the reduc- 
tion in Example 5.16 in the form of the section (Q/G) x (Q/G) G given by 
s(ro,ri) := p((ro, 0), ri). 
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In any case, the discrete unconstrained mechanical system associated to the 
reduced system by Theorem 10.7 is 

Q:=Q/G = M, with coordinate r, 

^d(ro,ri) = |((ri-ro)2 + (r?-r^)V4) 

with no forces since, by Example 5.6, the mixed curvature vanishes in the case 
b ~ 1, hence there are no reduced forces. Finally, the discrete Euler-Lagrange 
equation that determines the evolution of the {Q, Ld) system is 

(52) in+i - r,) - (r, - r,_i) + rk{{rl^, - rf) - {4 - rl^,))/2 = 0, 
that is, precisely, (22). 

Last, we adapt Theorem 7.1 to the reconstruction in the present setting. 

Theorem 10.10. Let G be a Chaplygin symmetry group of {Q.Ld^'D^'Dd). Let 
(go, 9i) G T^d o-nd r. be a discrete trajectory of the forced discrete mechanical system 
[Q/G, Ld,'D,'Dd, Fd) such that ■n{qj) = r.j for j = 0,1. Define the discrete curve 
q. in Q inductively by qq = qo and qk+i = hf{rk+i) for all k ^ N U {0} (here 
the horizontal lift is associated to the affine discrete connection whose horizontal 
space is Vd). Then q. is a trajectory of the original discrete mechanical system with 

qj = qj for j = 0, 1. 

Proof. The trajectory r. defines a trajectory {vk, rk+i) = yirkj^k+i) in G x (Q/G). 
Indeed, by Lemma 10.5 both curves satisfy simultaneously the corresponding vari- 
ational principles. Also, since qo £ 7r~^(ro), vq = p(qo,e) = p{qo, Ad{qa, qi))- In 
general, Vk = p{qk,e) for some qk eTT^^{rk). 

Then, since {qo, qi) € Vd, by Theorem 7.1, there is a trajectory q. in Q such that 
7r(gfc) = rk for all k and qo = qo, qi = qi- Furthermore, qk+i satisfies the formula 
in the statement because, qk+i satisfies (31) with Uk = e for all k. □ 

Example 10.11. We use Theorem 10.10 to reconstruct the trajectories of the 
system reduced in Example 10.9. Applying the recursive formula and (9) wc have 

qk+i = {xk+i,yk+i) = hf{rk+i) = {rk+i,yk + {rl+i - rl)/2), 

or, simplifying, {xk,yk) = {'''k,yo + (^^ ^ ^o)/^)' agreeing with the result obtained 
in Example 7.2. 

11. Reduced equations of motion: horizontal symmetries case 

In this section we specialize theorem 5.11 to the case where the original system 
is a mechanical system with horizontal symmetries. In the same way as in the 
Chaplygin case discussed in Section 10, in this case we go beyond the result of 
equivalence between the mechanical system in Q and a dynamical system in G x 
{Q/G) to obtain an equivalence between the discrete mechanical system in Q and 
another one in Q/G. 

Definition 11.1. Let M be a symmetry group of {Q, Ld, T), Vd). A closed subgroup 
G C A/ is said to be a horizontal symmetry subgroup for {Q, Ld, V, Vd) if 

(53) V^^q) n 2?g = V^{q) for ah qeQ. 
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From now on we will forget the group M and consider the action of G on the 
system. It is in this context that we specialize Theorem 5.11. Due to condition (53), 
we have that V*^ = S, the decomposition (2) of TQ becomes 

TQ = W® {0} ®S®H, 
u 

for any complementary subbundles % oi S m T) and W of 2? in TQ. Fixing one 
such decomposition we define a connection A on the principal bundle tt : Q — > Q /G 
requiring that Horj\^ = T-L. 

In the context of this section the discrete nonholonomic momentum map Jd 
defined by (32) has some special properties, which are studied next. 

Lemma 11.2. // G is horizontal symmetry subgroup the following statements are 
true. 

(1) 0^ = gx0 and {q^Y ^ Q X Q* . 

(2) Composing the nonholonomic momentum map Jd defined in (32) with the 
projection onto the second variable defines a momentum application Jd : 
Q X Q ^ Q* . Explicitly, for ^ € Q, 

Jd{(lo,qi)£. := -DiLd{qo,qi)^Q{qo)- 

(3) Any of the equivalent conditions of Theorem 8.5 is equivalent to the condi- 
tion that Jd is constant on the trajectory q. . 

Proof. Item 1 is a direct consequence of V"^ C while item 2 is clear from item 1. 

Last we check item 3. Let q. be a discrete curve in Q. Assume that equation (35) 
holds on q. for any section ^. Then, since any C S defines a (constant) section, 
evaluating (35) on this section yields [Jd{qk, qk+i) — Jd{qk~i , qk)) £, = 0. Thus, since 
^ G g is arbitrary, Jd is conserved on q. . 

Conversely, if Jd is constant on g. , equation (35) holds for constant sections of 
0^. But, it can be readily checked that if (35) holds for a section ^, it also holds 
for the section for arbitrary / : Q — > R. Then, since in our setup every section 
is a linear combination of constant sections with variable coefficients, we conclude 
that (35) holds for all sections. □ 

Below we construct an affine discrete connection adapted to the present geom- 
etry. Later we use that connection to specialize Theorem 5.11 to the horizontal 
setting. 

11.1. Affine discrete connection for horizontal symmetries. We recall two 
well known results. 

Lemma 11.3. If ^ e Q, g e G and q e Q, then CQilfiq)) = Zj«((Adg-i (e))Q(g)). 
Lemma 11.4. If g E G and qo,qi £ Q, then 

Jd{l'^'''^{qo,qi)) ^ Ad;-,{Jd{qo,qi)). 

When FiQxQ— >Risa smooth map, DiF : p\TQ — M is defined by 
DiF = dF\p'^TQ- Noticing that p*TQ ~TQ xQ, we have D2{DiF) : p*TQ R 
(where p2 ■ TQ x Q — > Q is the projection). As P2TQ ~ TQ x TQ, it is customary 
to consider D2DiLd : TQ x TQ — > R. It is easy to check that D2DiLd is bilinear 
in the tangent vectors. 
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Definition 11.5. Let id : Q x Q ^> M be a discrete lagrangian. We say that Ld is 
regular at (go, 91) G Q x Q if the bihnear mapping D2DiLd{qo,qi) '■ Tq^Q x Tq^Q — > 
E is nondcgencrate, i.e., if Xq € Tq^Q satisfies D2DiLd{qi)Tqi){XQ,Xi) — for aU 
Xi G Tq-^Q, then Xq — 0. In coordinates, the regularity condition becomes that 
the matrix ^ ^•'('io--<]i) ]-,e invertible. 

dqodqi 

Definition 11.6. Let G be a symmetry group of {Q, Ld, T>, Vd). We say that Ld is 
G-regular at (go, qi) Q x Q if the restriction of the bihnear form D2DiLd{qo, qi) : 
Tq^Q x TqjQ ^ M to V'^(go) x V'^(gi) is nondcgencrate. 

The notions of regularity introduced above have already been considered by other 
authors ([21], [10], [23]). In order to study the relationship between the regularities 
of a lagrangian and the fact that the discrete momentum values be regular we begin 
by recalling the following fact. 

Lemma 11.7. Let tt : Q — > Q /G a principal bundle and {yi, . . . , Vk} C g a linearly 
independent subset. Then, if q £ Q, {(wi)q((7), . . . , (wfc)Q(g)} C TqQ is a linearly 
independent subset. Furthermore, if the first set is a basis, then the second one is 
a basis too. 

Proposition 11.8. Let G be a horizontal symmetry group of [Q,Ld,'D,'Dd) with 
regular Ld, discrete momentum mapping Jd : Q x Q ^ q* and /i G g*. Then, 

• fi is a regular value of Jd and, consequently, :— Jj^^{fi) d Q x Q is a 
submanifold. 

• if for qo € Q we let Jf' '■ Q ^ Q* by J|°(gi) := Jd{qo,qi), is a regular 
value of and, consequently, {J'^)~^{P') C Q is a submanifold that, if 
not empty, has dimension dimQ — diniG. 

Proof. We have to prove that for all (go, 91) G Jfi, the map 

dJd{qo,qi) ■■ T(qa,qi){Q xQ) ^ T^0* ~ 0* 

is onto. Let {ei, . . . , e^} be a basis of g and {e^, . . . , e*} its dual basis. Hence 

-^d(go,'?i) = y] . (j)jiqo,qi)e*, 
where (j>j{qo,qi) = Jd{qo,qi){ej) = ~DiLd{qo,qi){ej)Q{qo). Then, 

dJd{qo,qi)iXo,Xi) = J2]^^iDl(l)j{qo,qi)iXo) + D2(l)jiqo,qi){Xi)) e* . 

Given i/j = ^j^'j G 0*, by the regularity of Ld there is Xi G Tq-^Q such that 

D2(l)j{qo,qi){Xi) ^ -D2DiLd{qo,qi){{ej)Q(qo),Xi) = aj 
for all j. Then, 

dJd{qo,qi){0,Xi) = 2_^^^^D2(l)j{qn,qi){Xi)e* ^ 2^^^^a.je* = V, 

so that dJd{qo, qi) is onto; thus /x is a regular value of Jd and standard results allow 
us to conclude that C Q x Q is a submanifold. 
We see in the previous computation that, for ?/' G g* 

dJ«°(gi)(Xi) = D2Jd{qo,qi)iXi) = ^, 
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SO that fi is also a regular value of Jf', hence, (J*) ^(/i) C Q is also a submanifold. 
The dimension of can be computed noticing that, being fi a regular value 

of J]°, if qi e iJj")-\^^), we have 

dim((J*)-i(/i)) =dim(ker(dJ]°(gi))) = dim(r,,Q) - dim(Im(dJj'' (gi))) 
= dim(g) -dim(G). 

□ 

Proposition 11.9. Let G be a horizontal symmetry group of {Q, Ldj'Dj'Dd) with 
Ld regular and G-regular. Then, for all qq ^ Q and fi G q* , if qi G 'g({(Zo}) H 
(J|")-1(a.), then T,,Q = T^J'^Hqo}) ®T,,{JfTHl-^)- 

Proof Assume that {J^y^il^) + 

since, otherwise, the statement is valid. As 
dim(Jj«)-i(/x) = dimg-dimG, it suffices to see that T,,/g({go})nT,,(J«")-i(^) ^ 
{0} for all qx as in the statement. If X-^ is in this last intersection, 

= dJf (<7i)(Ai) = -Z?2Z?iid(go,9i)((e,)Q(<Zo),Xi)e*, 
i. — 

so that D2DiLd{qa, qi){{ej)Q{qo), Xi) =0, for all j. By the G-regularity of Ld and 
being {(ei)Q(qo) • • • , (efc)Q(go)} a basis of TqJ^{{qo}), wc have that Xi^O. □ 

Definition 11.10. Let G be a horizontal symmetry group of (Q, Ld, T>, Vd) with Ld 
regular and G-regular. Given fi & Q* we say that G is a group of fi-good symmetries 
if, in addition, for each q G Q there is a unique g G G such that Jdi^'^il)) = In 
this case, we define ^ : Q ^ G hy ^{q) := .g. 

It is possible to extend the previous notion to systems where there are more 
than one g € G with the required property but, in this case, the action must be 
accompanied by a smooth unique determination of g. 

Proposition 11.11. Let G be a group of ^-good symmetries of (Q,Ld,'D,'Dd) for 
some /i G 0*. // /i satisfies Ad*{fi) = fi for all g G G, then <Z Q x Q defines an 
affine discrete connection in n : Q ^ Q/G of level 7, given by the fi-goodness of G. 

Proof. By Proposition 11.8, is a submanifold of Q x Q and, by Lemma 11.4, 
it is also G-invariant. According to the definition of 7, Jd{q,l'^i^g-^{q)) = M; so 
that r C J7^. Last, by Proposition 11.8, J^i{q) C Q is a submanifold and, using 
Proposition 11.9, we conclude that defines an affine discrete connection. Notice 
that 7 is G-equivariant by Lemma 11.4 and the condition Ad*{^) ~ /i for all 
gGG. □ 

Remark 11.12. In the context of Proposition 11.11, the condition that = ji 

for all g £ G can be avoided by considering the (probably smaller) symmetry group 
{ff S G : Ad*{^) = fi} instead of G. Indeed, G^ is a Lie group and a 
symmetry group of {Q, LdjV^Vd)- Since G satisfies (53), V^iq) C Vq for aU q, so 
that V^f^iq) C Vg for all q. Hence, (53) holds if wc put G^ instead of M and G 
and run the reduction arguments in this setting. This same remark applies to the 
statements made in what remains of this section. 
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11.2. Reduction. In this section we assume that G is a group of /i-good symme- 
tries of a regular and G-regular system for some ^ € q*. In the previous section 
we saw that if /i satisfies Ad*{fi) ~ fi for ah g G G, we can define an affine discrete 
connection Ad such that Horj\^^ — J^. Below, we use this connection to specialize 
Theorem 5.11 to the context of horizontal symmetries. 

As in the Chaplygin case studied in Section 10.1, using y we define a forced 
discrete mechanical system on Q/G with discrete lagrangian Ld(ro, ri) := y*{Ld) = 
Ld ° y, variational constraints T) := dTr{'D), kinematic constraints Vd ■= y^^iT^d) 
and forces F ■.= y*{F). 

Theorem 11.13. Let G be a fi-good symmetry group of (Q, Ld^'D^'Dd) for some 
fi E Q* . Assume that Ad*{fi) = fi for all g G G. Let q. be a discrete curve in Q and 
Tfc := 7r((7fc) the corresponding discrete curve in Q/G. Then, if JdiQk, Qk+i) = fJ- for 
all k, the following statements are equivalent. 

(1) {qk, qk+i) £ T^d for all k and q. satisfies the variational principle dSd{q){Sq.) = 
for all vanishing end points variation Sq. such that dqk S "Dq^ for all k. 

(2) {qi;,qk+i) G T^d for all k and q. satisfies the discrete Lagrange-DAlembert 
equations (3) for all k. 

(3) r. satisfies the variational principle 

dSdir.)i6r.) -^^J^'(F-(rfe,rfe+i) + F+(rfe_i,r,))(5r,) 

for all vanishing end points variations such that 6rk G "Dr^ for all k. In 
addition, (rk,rk+i) G Vd for all k. 

(4) r. satisfies (51). In addition, (rk,rk+i) G "Dd for all k. 

Proof. In the present context Jd is conserved over q.. Hence, equation (35) is 
satisfied and, by Theorem 8.5 the vertical equations in Theorem 5.11 are satisfied. 
Then, by Theorem 5.11, points 1 or 2 of the present result are equivalent to either 
one of 

(i) Item 3 in Theorem 5.11 holds for all horizontal variations, (^^ = for all 
k). 

(ii) For aU k {vk,rk+i} G Vd and condition (16) holds. 

Notice that Wk = Ad{qk, qk+i) ~ e due to the conservation of Jd and the choice 
of Ad- Then, (ufc,''A;+i) = (piqk, e),rk+i) = 3^(rfe,rfe+i) G t>d if and only if 
{rk,rk+i) G y~^{T>d) = "Dd- Thus the constraint conditions that appear in items i 
or ii are equivalent to the ones that appear in 3 or 4. 

By Lemma 10.5, the variational principle that appears in item i is equivalent to 
the variational principle in point 3. Similarly, by Lemma 10.6, the equations that 
appear in item ii are equivalent to the ones that appear in 4. □ 

The following reconstruction result is the analogue of Theorem 7.1 which is valid 
in the current setting. 

Theorem 11.14. Let G be a horizontal symmetry group of {Q,Ld,T>,T>d) and 
{qo, qi) G Vd. Let /i ~ Jdiqo, qi) G 0*- Assume that Ad*{^) = /i for all g G G and 
that G is ^-good. 

Let r. be a discrete trajectory of {Q /G,Ld, T>, Vd, Fd) such that ^{qj ) = rj for j = 
0, 1. Define the discrete curve q. in Q inductively by qo = qo, arid qk+i = h'^diTk+i) 
for all k (here the horizontal lift is associated to the affine discrete connection whose 
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horizontal space is J^l). Then q. is a trajectory of the original discrete mechanical 
system that satisfies qo = qa and qi = qi. 

Proof. See the proof of Theorem 10.10. □ 

11.3. An example of reduction of horizontal symmetries. In this section 
we apply the analysis developed in the previous sections to a system that exhibits 
horizontal symmetries. The system is a discretization of the nonholonomic free 
particle considered by J. Cortes on page 100 of [11]. More precisely, the system 
{Q,Ld,'D,'Dd) has Q := and 

L4<10,qi) f - q^o f + (^f - Q^o)' + ill - lo f), 
where q := {q^ ,qy ,q^). The constraints are 

:= {xd^+ydy + zd, eT,Q:y = q^i} = {{d,\^ + q^dy\^),d,\^), 

Vd := {(go,gi) eQxQ-.ql-q^i^ {{qtf - {q'^f)/2]. 

The group M := R? acts on Q by l'^{q) :~ {q'^^q^ + 5^,5^ + 5^), where g 
(5^1 ff^) S M. The Lie algebra of M is m that we identify with M^. The cor- 
responding lifted action is l^^{q,v) = {l^{q),v). Therefore, Ld, V and Vd are 
Af- invariant. 

The vertical space for the action is V^' = (^^l,' ^-Ig) '^'^qQ^ ^ti^'t 

Thus, the closed subgroup G := {0} x R C M, satisfies Vg n V^' = Sg = , for 
all q £ Q. Hence, G is a horizontal symmetry (sub)group of the system. We have 
g = {0} X R c ^ m. 

Clearly, Ld is regular as well as G- regular. Also, being G abclian, Ad*{fi) — ^ for 
all g € G and /i G g*. A simple computation shows that Jd{qo, qi) = ra{ql — (7o)l*, 
where 1* denotes the basis of g* that is dual to (0, 1). Then, if g G Q and ^ = /i^l*, 
J'di^'giq)) — M has a unique solution g := (0, ^M^) G G and G is a group of fi- 
good symmetries. Hence Theorem 11.13 applies to identify the reduced system 
with a forced discrete mechanical system on Q/G ~ M^. Below we give an explicit 
description of this reduced system. 

As in the previous sections, we use the affine discrete connection Ad whose 
horizontal space is, for a fixed /i G g*, Jfj, = {{qo, qi) E Q x Q : m{qf — q^) = ^^}. 
Equivalently, Adiqa^qi) = (0,qf ^ 9o ^ m^^) ^ '^^'^ discrete horizontal lift of 
^rfis7^f(r) = (r',r",i/.- + gg). 

As an intermediate step we have to describe the reduced system on G x [Q/G). 
Since 7r:Q— s^Q/Gisa trivial principal bundle with structure group G we apply 
the description of the reduced system given in Section 9. The resulting system is 

'G x (Q/G) ~ (Q/G) X G X (Q/G) = R^ x R x R^, 
^ L* (r,,^fe, r,+i) = f ((r;,+i - r',f + {r'l^, r'[f + {^f,^ + duf), 

Vl = {dr'\r+r'dr''\r) C TriQ/G), 

V\ = {(r,,^fc,r,+i) G (Q/G) x G x (Q/G) : = {{r'^^.f K?)/^- 
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In order to complete the description of the reduced system we need to compute the 
reduced forces Fd- We fix a sphtting 




= Wg =U, =5, =«, 



which determines the nonholonomic connection A. In particular, notice that since 
dAd{qo,qi) = (0, d^^lqi — dz|<jo) ^-nd Horj\{6q) G {dx\^ + q^dy\^) the mixed curvature 
Bm vanishes, and so docs the reduced force Fd- 

The inclusion y is y(rk,rk+i) {r^, e,rk+i) so that the reduced mechanical 
system {Q, Ld,'D,'I)d) has Q — Q/G and 

(54) }'Dr = {dr'\^ + r'dr"l), 

[Vd = {(rfe,r,+i) G (Q/G) x (Q/G) : r'^^, - <f = (iki);^!!!)!}. 
The discrete Lagrange-D'Alembert equations (3) for this system are 

-K+i-^;) + K.-'i.-i) = A.r;, 

rUi-^k=iirUiy-i^k)')/^- 
Those equations can be easily solved to determine the evolution of the reduced 
system, r.. The corresponding trajectory q. of the original system is obtained 
recursively, according to Theorem 11.14, by, for a given qk, defining 

qk+i Hqk+i,ql+i,qk+i) = hf{rk+i) = iri+l,4+l,^^i'' + qk) 
={r'k+i^r'l_^i,ql~qo+ql). 

or, 

{ql. ql ql) = K.rl k{ql ~q^) + q^). 
Interestingly, the discrete mechanical system {Q^Ld^V^Vd) defined by (54) still 
has a residual symmetry group. The group M/G ~ M acts on Q/G via the action 
induced by l'^ on Q/G. This is precisely the system and symmetry group whose 
reduction and reconstruction as a Chaplygin system was discussed in Examples 10.9 
and 10.11. Using those results we get the trajectories of the original system 

{<ll^<ll<ll) - {rk, {rl rl)/2 + j/o, k{ql q^,) + q^), 
where is determined from (52) with tq = q^ and ri ~ qf. 

Remark 11.15. The presence of this second reduction step by M/G is an example 
of reduction by stages for discrete mechanical systems, a topic that will be discussed 
elsewhere. 
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